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ABSTRACT 



Besides the conventional hadrons containing valence quarks and valence an- 
tiquarks, quantum chro mo dynamics (QCD) suggests the existence of the hybrid 
hadrons containing valence gluons in addition to the quarks and antiquarks, and 
some experiments may have found some. A decisive experimental confirmation of 
its existence, however, is still needed. At present, lattice simulations have offered 
the practicable ways of theoretically guiding us to search for the hybrid states. 

In this dissertation, we study the spectroscopy and the decay rate of the heavy 
hybrid mesons with a heavy b quark, a heavy b antiquark and a gluon (bbg) to the 
chosen channels, and use lattice methods to calculate the transition matrix elements 
in full QCD. We are particular interested in the spin-exotic hybrid mesons. 

For sufficiently heavy quarks (e.g., b quark), we use the leading Born-Oppenheimer 
(LBO) approximation to calculate the static potential energy at all bb separations. 
Then, by solving the Schrodinger equation with this potential, we reconstruct the 
motion of the heavy quarks. In a similar way we can determine decay rates. In 
this dissertation, we use the lattice method to calculate the mass of the /o meson 
at a single lattice spacing and hght quark mass, namely, m/g = (768 ± 136) MeV. 
Most of all we consider the decay channels involving the production of a scalar 
meson. We obtain the partial decay rate (F) for the channel H — >■ Xb + + t^, 
namely, F = 3.62(98) MeV. All of our results are consistent with those of other re- 
searchers. Knowledge of the masses and the decay rates should help us considerably 
in experimental searches for the hybrid mesons. 
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CHAPTER 1 
INTRODUCTION 

Besides the conventional hadrons, quantum chromodynamics (QCD) predicts 
the existence of the hybrid hadrons, which are hadrons but also contain the excited 
gluon fields. Early results in e+e" annihilation cross section from the CUSB and 
CLEO collaborations [1, 2] showed that there exists a complex resonance structure 
in the energy region above the T(45') (near 10.7 GeV), where the possible lowest 
hybrid states (bbg) are predicted in Ref. [3]. However, because it is difficult to 
differentiate a resonance from a threshold effect, we just have a "hint" of the 
existence of the hybrid states, a decisive experimental confirmation is still needed. 

We need some theoretical guidance to search efficiently for the hybrid states. 
Due to the nonperturbative nature of the gluons, however, a perturbative theory 
in this field is still not available. Hence, a clear understanding the hybrid states by 
analytical methods is very difficult. At present, the lattice simulations have offered 
practical ways of theoretically guiding us to search for the hybrid states. 

The heavy hybrid mesons can be investigated directly by the leading Born- 
Oppenheimer (LBO) approximation. The results from our LBO approximation can 
compare with the corresponding results from other people's experiments. We focus 
on the heavy hybrid mesons (QQg), which consist of a heavy quark (Q) and a heavy 
antiquark (Q) plus a gluon (g). We study its spectroscopy and the decay amplitude 
from the heavy hybrid meson to {QQ){qq), and use lattice methods to calculate the 
transition matrix elements in full QCD [4]. Finally, we study its phenomenological 
implications. We are particular interested in the spin-exotic hybrid meson, which is 
the meson with J^*^ quantum numbers which are not permitted in the current quark 
model. The lightest one is beheved to be with quantum numbers J-^^ = [5], 



2 



which does not mingle with a nonhybrid QQ bound state [4]. 

The spectroscopy of the hybrid mesons (including spin-exotic hybrid mesons) 
has been extensively studied over the past years [5, 6]. The BO approximation was 
proposed for the study of the hybrids states in Refs. [7, 8]. For sufficiently heavy 
quarks (e.g., h quark), the LBO approximation justifies a treatment in which the 
heavy quark {Q) and antiquark {Q) are regarded as approximately static, and we 
assume that the gluon field adjusts rapidly to their movement. Hence, for various 
gluonic excitations of the flux tube of QQ, we compute its corresponding static 
potential energy at QQ separation R. Then by numerically solving the Schrodinger 
equation with this potential, we reproduces the motion of the heavy quarks. In a 
similar way we can calculate the decay rate. 

Decay modes are also intensively investigated in Rcf. [4]. However, at present, 
our knowledge about the decay channels is still limited. We make a contribution to 
the understanding of the spectroscopy of the hybrid mesons and its decay model. 
There are many new features in our simulation: (1) We adopt gauge configurations 
with improved actions, which include the Symanzik improvements at 1-loop level, 
and "Asqtad" fermion action [9]. (2) We choose the improved wave functions for 
the spin-exotic gluon channels, which are used to integrate all the decay rate at all 
QQ separation R. (3) We consider decay channels — )■ Xfc + ^ + ^ involving the 
emission of two pions (here H stands for hybrid exotic state) . (4) The mass of the 
dynamical sea quark we choose is closer to the physical value. 

Our project needs three classes of adjustable input parameters: 1) lattice spac- 
ing: a, 2) light quark masses {u,d,s quarks), 3) the mass of the bottom quark 
(6 quark). These input parameters are obtained from experiment. We obtain the 
masses of u, d, s quarks from the masses of n, K mesons, and we obtain the lattice 
spacing a from the splittings of T(25') — T(IS') and T(1P) — T(IS') in T spectrum. 
And from the B meson, we obtain the mass of the bottom quark [10]. This project 
help us gain the knowledge about "soft gluons." Furthermore, we can predict the 
mass of some mesons, e.g., the /o meson. Most importantly, we can calculate the 
partial transition rate F, which has been poorly understood. Knowledge of the 
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masses and decay rates helps us in experimental searches for the hybrid mesons. 

In this dissertation, we use the MILC version of the Kogut-Susskind staggered 
quark action, and the gauge action is the Liischer Symanzik improved action. 
Also we use gauge configurations, which were generated by the hybrid-molecular 
dynamics "R algorithm" [11]. These configurations contain three flavors: two light 
quark {u, d), and a strange (s) quark. In our simulation, the light quarks are 
degenerate in mass, and we choose the bare strange quark mass and up, down 
quarks mass to be close to the corresponding physical values (i.e., arriu — 0.005, 
and anis — 0.05). Prom Ref. [12] we find the physical mass of the pion meson is 
about 140 MeV, from our lattice calculation, we get the mass of the pion meson is 
about 230 MeV. The lattice space a is in our simulation of around 0.12 fm [10]. 

This dissertation is organized as follows: In Chapter 2, wc provide background 
material on particle physics. We first give a brief review of the quark model. 
Then we briefly discuss a few conservation laws. Finally, we discuss about the 
quarkonium (heavy quark, antiquark states). In Chapter 3, we provide background 
material about QCD. First we give an overview of continuum QCD. Then we give 
an overview of lattice QCD. In Chapter 4, we use the LBO approximation to study 
the spectroscopy of the hybrid mesons. We first study the method of the LBO 
approximation. Then we introduce the static potential. Finally, we calculate hybrid 
quarkonium states through the LBO approximation. In Chapter 5, we first briefly 
review the hybrid states on the lattice. Then we give an overview of hybrid meson 
decays. Finally, we give a overview of the decays from the lattice. In Chapter 6, 
we explain our calculational method in detail. In Chapter 7, we first use time- 
dependent perturbation theory to study the Chab correlator. Then we introduce 
a method to extract transition amplitude x directly from lattice. In Chapter 8, we 
give all the lattice numerical simulation results. First we present the calculation of 
the mass of the /o meson. Then we show the detailed procedure to obtain lattice 
transition amplitude x, and calculate the decay rate. In Chapter 9, we first discuss 
the evidence for the /o meson and the hybrid exotic mesons. Then we discuss hybrid 
exotic quarkonium decay. Finally, we summarize our results. 



CHAPTER 2 
INTRODUCTION TO PARTICLE PHYSICS 

Particle physics is the study of the fundamental particles that make up matter, 
and how they interact with each other [13]. At present, the "quarks," along with 
"leptons," are believed to be fundamental, and the quark model is a popular theory 
in this field. 

In this chapter, we provide background material about particle physics. We first 
give a brief review of the quark model. Then briefly discuss a few conservation laws. 
Finally we discuss about quarkonium (i.e., the bound states of heavy quarks and 
antiquarks). We focus on just a few details needed for this dissertation. Several 
good textbooks [13] [14] provide comprehensive background material. 

2.1 Quark model 

In particle physics, we define a hadron as a subatomic particle which takes part 
in the strong nuclear force. In the current quark model, hadrons are made up of 
quarks, antiquarks, and gluons. Usually hadrons are further classified into many 
types (e.g., baryons and mesons, etc). In this dissertation we are interested only in 
baryons and mesons. 

In 1961 Murray Gell-Mann [13] introduced the Eightfold way to arrange the 
hadrons which were discovered at that time. The eight lightest baryons (n, p, E+, 
E°, A, S°) and the eight lightest mesons X+, , K^, 7r°, 7r°, 7r+, 
rf) fill a hexagonal array, with two particles at the center, according to their charge 
and strangeness. They are called a baryon octet, and meson octet, respectively. ^ 



^See the figures on page 33, 34 in Ref. [13]. 
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In order to explain and understand the Eightfold way, in 1964, Murray Gell- 
Mann and George Zweig independently proposed [13] that just three fundamental 
constituents (and their antiparticles) combined in different ways could explain the 
Eightfold way. Murray Gell-Mann called these elementary constituents "quarks," 
and the three types (or "flavors") of quarks were named up (m), down (d), and 
strange (s), respectively [13]. 

In the late 1960s the evidence for the existence of "quarks" became clear. In 
1974 a new particle J/^ [13] was discovered, which was later understood to be a 
bound state of a new quark-antiquark pair (cc). The new fourth quark is known as 
charm (c). In 1977 another new meson called the Upsilon (T) [13] was discovered at 
Fermilab, which was later shown to be the bound state of the new quark-antiquark 
pair (bb). The new fifth quark was called bottom (6). Finally, in 1995, the evidence 
for the existence of the "top" (t) quark was obtained at Fermilab. Hence, the 
current quark model with six "flavors" was fulflUed. 

In the current quark model, we have the six quarks: up, down, strange, charm, 
bottom and top. By convention we usually label them by their flrst letters: u, d, s, 
c, b and t. Moreover, each quarks has its partner antiquarks. We label them u, d, s, 
c, b, t, respectively. The charges of up, down, strange, charm, bottom (or beauty), 
and top (or truth) quarks are |e, — |e, — |e, |e, — |e, |e respectively, where e is 
the charge of the electron. For each quark (g), the corresponding antiquark (g) has 
the opposite charge. Hence, the charges of antiup {u), antidown (d), antistrange 
(s), anticharm (c), antibottom (6), antitop (t) quarks are — |e, |e, |e, — |e, |e, 
— |e, respectively [13]. 

The six types of quarks (phis their six antiquarks) are classifled into three pairs, 
which are listed in Table 2.1. They are the up-down, the charm-strange, and the 
top-bottom multiplets. The top block is the first generation, the middle block is 
the second generation, and the bottom block is the third generation [13]. 

Like quarks, there are six types of the leptons [13], which also are classifled into 
three pairs: electron-neutrino (e — i/g), muon-neutrino {/i — i/fj), and tau- neutrino 
(r — Ut-). The electron, muon, and tau each carry a negative charge, but the three 
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Table 2.1. Three quark generation. Where Q is the electric charges, D is the 
"downness," U is the "upness," S is the "strangeness," C is the "charm," B is the 
"beauty," and T is the "truth." 
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neutrinos (z/g, ^t) carry no charge. 

The quark model assumes [13] that the quarks bind with each other to form 
the particles. There are two basic types of the combination. The first type of 
the combination is "baryons," which are made up of three quarks. (Similarly, an 
antibaryon is made up of three antiquarks). Listing all the combination of three 
quarks gives the baryon decuplet [13]. One example of a baryon is the proton (p), 
the proton is made up of "uud." ^ whose charge is +1. 

Another type is the meson, which is made up of a quark and an antiquark. 
Counting all the quark- antiquark combinations yields the mesons in which we are 
especially interested. We list these mesons in Table 2.2, the top block is the meson 
nonet, the middle block contains "charmed mesons" [12], and the bottom block 
contains "bottom mesons" [12]. The "charmed mesons" and "bottom mesons" are 
heavy. 

The quarks and antiquarks carry spin |. When a quark and an antiquark bind 
to form the mesons in a state of zero orbital angular momentum (L = 0), there 
are spin-0 combinations that are called "pseudoscalar meson" (i.e., tt"*", tt", 7r°, K'^, 
K^, K~ , K^, 77, 77'), and spin-1 combinations that are called "vector meson" (i.e., 
p+, p~, K*^, K*^, K*~, (f), u). In the language of group theory, each 

quark can occur as spin-up (mc,. = |) or spin-down (m^ = — |) states, and two spin 
states belong to the fundamental representation of spin SU{2) (2). Similarly, each 



^Here we just include the valence quark. Later in this chapter we will consider the sea quark, 
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Table 2.2. The top block is the meson nonet, the middle block is "charmed 
mesons," and the bottom block is "bottom mesons." Here we just include the 
mesons which we are interest in this dissertation. Where the Q stand for charge, 
and S is strangeness, the C stand for charm, and B is bottom (or beauty). 
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antiquark can also occur as spin- up (mg = |) or spin-down (m^ = — |) states, and 
two spin states belong to the complex conjugate fundamental representation of spin 
SU{2) (2). We combine two two-dimensional fundamental representations (2, and 
2) of spin SU{2) to obtain a spin triplet (3, spin s = 1) and a spin singlet (1, spin 
s = 0), 

20 2 = 30 1. (2.1) 

Hence, the spin singlet state is pseudoscalar mesons (namely, the total spin s = 0), 
and the spin triplet state is vector mesons (namely, the total spin s = 1) [13]. 

In addition to the electric charge, quarks also carry "color charge" [13]. In the 
language of group theory, each quark takes "red"(r), "green" (gr), or "blue" (6), and 
the three color states belong to the fundamental representation of color SU{?>) (3). 
Similarly, each antiquark occurs as "antired" (f ) , "antigreen" (^), or "antiblue" (6) , 
and three anticolor states belong to the complex conjugate fundamental repre- 
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sentation of color SU{3) (3). We combine two three-dimensional fundamental 
representations (3, and 3) of color SU{3) to form a color octet (8) and a color 
singlet (1), 

3(8)3 = 80 1. (2.2) 

The quark model asserts that all macroscopically isolated particles must be "color 
singlet." For example, mesons are the combinations of a quark of color and a quark 
of its anticolor to form a neutral color charge. 

The force between color charged particles is the "strong force." The carriers of 
the strong force are gluons [13], which carry no electrical charge, but have both a 
color and an anticolor charge. According to Eq. (2.2), the gluon can exist as color 
octet states or color singlet state. However if the gluon of color SU (3) singlet state 
exists as mediator, it can also occur as a free particle. An interesting fact about 
quarks and gluons is that we have never been observed a macroscopically isolated 
quark or gluon. For this reason, the color-charged quarks (or gluons) are believed 
to be confined within a hadron [13]. Hence, we have eight physical gluon states. 

By convention valence quarks for a given hadron arc the quarks (or antiquarks) 
that give rise to the quantum numbers of the given hadron. For example: the 
proton's quantum numbers are characterized by two "m" and one "d" quark valence 
quark. At any given momentum the proton may contains extra uu pairs, or dd pairs, 
or even ss pairs, cc pairs, etc. We call these extra quarks sea quarks [13]. 

2.2 Some conservation laws 

The study of the interactions and the decays has resulted in a few conservation 
laws. For the strong interactions and decay, the conserved quantities include baryon 
number (A), lepton number (L), upncss (U), downncss (D), strangeness (S), charm 
(C), bottom {B), top (T), isospin (/), parity (P), and charge-conjugation (C). 
These conservation laws, together with the classical conservation laws (such as the 
conservation of energy, momentum, charge) should apply in the strong interactions 
and decays. 
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In this section, wc just focus on the conservation laws of parity (P), isospin (/), 
and charge-conjugation (C), which are needed for this dissertation. 

2.2.1 Parity 

Many physical processes like strong interactions have a property of parity invari- 
ance (or "mirror symmetry"). This means that the mirror image of any physical 
process can also ocur with the same probability. In this subsection, we repeat some 
discussion in Sec. 4.6 in Ref. [13]. 

In physics, the parity basic operator is an inversion operation, in which the 
object is flipped to the opposite location through the origin [13]. 

Quarks have an intrinsic parity, which is defined to be +1, and for an antiquark 
the parity is —1 [13]. Since the quantum number parity for the composite system 
is multiplicative, the intrinsic parity of the baryon octet and decuplet is +1, and 
the intrinsic parity of the meson is —1 [13]. For the excited states of the mesons, 
the parity is given by [13] 

p = {-iy+\ (2.3) 

where I is the orbital angular momentum of the meson. 

The parity is not conserved by the weak force, with the result that all the 
neutrinos are found to be "left-handed" [13]. 

2.2.2 Isospin 

Isospin is a physical term that was introduced to describe a group of particles 
that has almost the same mass [13]. For example, the doublet of the proton and 
the neutron is said to have isospin ^ {I — |), with the third component /s = +| 
for the proton and — | for the neutron. Another example is the three pions (7r"'",7r'^, 
7r~), which compose a triplet. Hence, its isospin is 1. The third components are 
-1-1 for the tt"*", for 7r°, and —1 for tt", respectively. 

In 1932 Heisenberg suggested that isospin is conserved in strong interactions. In 
the language of group theory, we can declare that strong interactions are invariant 
under an internal SU{2) symmetry (or internal isospin symmetry) [13]. Hence, 
nucleon (p, and n) is the two-dimensional representation of the SU{2) (/ = i), and 
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pion meson (tt"*", 7r°, and 7r~) is the three-dimensional representation of the SU{2) 

Isospin for the particles is related to other quantum numbers for the particles 
by general Gell-Mann-Nishijima formula [13] 

^ A + U + D + S + C + B + T 

Q = ^ , (2-4) 

where Q is the charge of the particle, A is baryon number, U is upness, D is 
downness, S is stangeness, C is charm, B is bottom, and T is the top [14] . 

At the current quark model, the bare masses of the up and down quarks are close. 
Hence, the "up" and "down" quarks can form an isospin doublet (/ = |), and all the 

other flavors {s,c,b,t) carry isospin zero [13]. By convention the third component 
of isospin of the up{u) quark is assigned to be |, and the third component of isospin 
of the down {d) quark is assigned to be — | (like proton and neutron) [13]: 



u 



2 2/ 



The antiup {u) and antidown {d) quarks also form an isospin doublet [13]: 



u 



2 2 



5 5>- P-^) 



By convention the third component of the isospin of u quark is assigned to be — | 
and the third component of the isospin of the d quark is assigned to be | [13]. 
If we add two particles whose isospins are / = |, we obtain an "isotriplet" [13] 

|1 1) ^ -ud, |1 0) = ^^-^"^ ^ \l-l) = du, (2.7) 

V2 

and an "isosinglet" 

uu + dd , , 

|0 0) = — (2.8) 

In the language of group theory, the u, d quarks belong to the fundamental 
representation of isospin SU{2) (2), and the u,d quarks belong to the complex 
conjugate fundamental representation of spin SU{2) (2). We combine two two- 
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dimensional representations (2, and 2) of spin SU (2) to obtain an "isotriplet" (3, 
isospin / = 1) and an "isosinglet" (1, isospin / = 0), 

2®2 = 3©1. (2.9) 

For the peudoscalar mesons the isotriplet is the pion. For the vector mesons 
the isotriplet is the p meson. Hence, the quark content (or the wave function of 
"flavor") for 7r°(or p°) meson is given by [13] 

A P° = (2.10) 
Isospin is connected with a conservation law, as shown by the strong process 



p + p > d + TT"*" 

I - - 1 (2.11) 

2 2 ^ ^ 

Here we assume the isospin of the deuteron (d) is consistent with the observation 
that there are no other bound state of two nucleons [14]. It is obvious that the 
above process conserves charge, angular momentum, baryon number, and isospin. 

2.2.3 Chcirge conjugation 

Charge conjugation is a physical operation that replaces each particle in a process 
with its the corresponding antiparticle [13] . It changes the sign of the charge, baryon 
number, lepton number, strangeness number, charm number, beauty number, and 
truth number of the particle. It does not, however, change the mass, energy, 
momentum, and spin of the particle [13]. 

We can show that for a system of a quark (spin-| particle) and its antiquark, 
the charge conjugation (C) is given by [14] 

C = (-l)'+^ (2.12) 

where / is total orbital angular momentum of the meson, and s is total spin of the 
meson. 

If we consider the ground states (/ = 0) of mesons, then for the pesudo-scalar 
mesons {s — 0), hence C = +1; for the vector mesons {s — 1), thus C — —1. 
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One kind of the notation for these states is to indicate their total angular 
momentum {j — I + s), parity, and charge conjugation explicitly. For the ground 
states {I — 0) of the pseudoscalar meson, the total angular momentum is zero (i.e., 
j — + — 0), and they have negative parity and positive charge conjugation (C). 
Hence 

jPC ^ Q-+ ^2.13) 

For the ground states of the vector mesons, s = 1, we obtain j — 1, and they have 
negative parity and negative charge conjugation (C). Hence 

J^^ = 1—. (2.14) 

Charge conjugation is conserved for strong and electromagnetic interactions, but 
not for weak interactions. We know that charge conjugation turns a left-handed 
neutrino into a left-handed antineutrino, and changes right-handed neutrino into a 
right-handed antineutrino. In weak processes, there exist only left-handed neutrinos 

and only right-handed antiucutrinos. Hence, it is obvious that charge conjugation 
is not conserved in weak processes [13, 14]. 

2.3 Quarkonium 

In the quark model, we regard all the mesons (51^2) as two-particle bound 
states [13]. The potential model works well for hydrogen and positronium. It 
is reasonable to ask whether the potential model works for the mesons or not [13]? 

After the discovery of the bound state cc (or charmonium states), potential 
models of QCD have become an important tool in understanding quarkonium (i.e., 
the bound state of a heavy quark and heavy antiquark, e.g., T{bb), cb, etc.) [13]. 

In hydrogen and positronium model, there exists only Coulombic potential [13]. 
(In QED, the force between electrical charges is mediated by the exchange of 
photons.) However, for the meson, we know that the force between color charged 
particles through gluon exchange is very powerful. Thus, the short distance be- 
havior of QCD is dominated by this "strong force." Hence, the potential of a 
meson should be Coulombic potential at short range [13]. However, due to quark 
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confinement, the potential can increase indefinitely when the distance r between 
quark and antiquark becomes large. Here we assume no qq pair production. At 
this range, it is standard to choose a linear fit (i.e., V{r) ~ r) (and it is confirmed 
in numerical simulation). Thus, in this dissertation, we choose the overall potential 
to be [13, 14]. In Chapter 4, we will note that for the potential, we choose a 
different model. 

V{r) = -— + k2r, (2.15) 
r 

where ki, k2 are just two constants. 

The mesons, which are made up of light quarks (i.e., u,d,s), are intrinsically 
relativistic [13]. The mesons, which are made up of heavy quarks (i.e., b,t), are in- 
trinsically nonrelativistic [13]. The classical quantum mechanics can describe these 
systems pretty well. From Fig. 5.7 in [13], we can see that potential model describe 
the charmonium very well. For the bound state bb (also called bottomonium) , the 
potential model describes it very well also. Please see Fig. 5.9 in [13]. 

In this dissertation we study the hybrid quarkonium (or hybrid meson) . Hybrid 
mesons are the bound state of QQ with nontrivial excited gluonic components 
(QQg), where Q is a heavy quark (e.g., b, t quark), and Q is a heavy antiquark. 
A hybrid spin-exotic state is a hybrid meson whose J^*" value is permitted in the 
quark model [4]. From the current quark model, the quantum numbers for QQ 
composites are given in Table 2.3. 

Because the quantum numbers J^^ = 1"'*', 0"^" and 2+~ cannot exist for the 
meson in the quark model, they are hybrid spin-exotic mesons. In Chapter 4, we 
give more details about hybrid spin-exotic mesons and the potential model. 



Table 2.3. The ^^+'^Lj or J^^ for QQ. 



L 


Singlet (5 = 0) 


Triplet (5 = 1) 





^So(0-+) 


%(o--) 


1 


ip,(l+-) 


%i,2(0++,l++,2++) 


2 


^D2(2-+) 


3Di,2,3(l--,2--,3--) 


3 


iF3(3+i) 


=^F2,3,4(2++,3++,4++) 



CHAPTER 3 



LATTICE QCD 

It is well known that quantum electrodynamics (QED) describes the interaction 
of the charged particles and photons. In the same way, quantum chromo dynamics 
(QCD) is the popular theory of the strong interactions of the colored particles, which 
is formulated in terms of the quarks and gluons. QCD has been very successful in 
predicting the phenomena of high energy hadron physics. 

Quantum chromodynamics is a well-developed theory, but calculations are very 
difficult, particularly at energies below about 1 GeV, where the strong coupling 
constant is of the order of unit, so perturbation methods fail. Fortunately, lattice 
QCD provides a nonperturbative method in this field. Lattice QCD is a quantum 
field theory that formulates QCD on a discrete Euchdean spacetime grid. We 
explain later. The main benefits are that many calculations in QCD take on forms 
famihar in classical statistical mechanics. Since lattice QCD does not introduce 
new field variables, it inherits the basic characters of QCD. 

In this chapter we provide background material about QCD. We first give an 
overview of continuum QCD. Then we give an overview of lattice QCD and some 
special ways to incorporate fermion fields. Finally, we explain how to generate 
the configurations. We focus on some details needed for this dissertation. Several 
textbooks [15, 16] provide comprehensive background material. 

3.1 Overview of continuum QCD 

In this section, we first introduce the full gauge invariant classical Minkowski 
Lagrangian for QCD by directly starting from the Minkowski Lagrangian of a free 
Dirac field for a given fiavor, and then introducing the gluon field through the 
help of the parallel transport. Then we write the path integral representation for 
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the physical vacuum expectation values of any observables. Later we change this 
path integral representation to the corresponding path integral in the Euclidean 
space-time. Also we associate the path integral formulation of the partition function 
in QCD with the corresponding partition function of quantum statistical mechanics. 

3.1.1 Continuum gauge field theory in QCD 

On the lowest level, QCD describes the interactions of the quarks and gluons. 
In Chapter 1, we explained that the quarks {u, d, s, c, b, t) and the antiquarks {u, 
d, s, c, b, t) are fermions of fractional electric charge and color charge. Each of the 
quarks carries a SU (3) color index. According to the current colored quark model, 
the three colors of a given flavor have the same mass, although different flavors 
have different masses [13]. Hence, the wave function ipc of each quark contains a 
triplet of fields {c — 1, 2, 3). The wave function t/j^ of each antiquark also contains 
a triplet of fields [16] [13]. 

Now we introduce the gauge-invariant Lagrangian in continuum QCD in the 
Minkowski space. The starting point is the Minkowski Lorentz invariant Lagrangian 
of the free Dirac field for a given fiavor [13] [15] 

3 

= E^c(a;)(^7'^a^ - m)M^), (3-1) 

c=l 

where c is the quark color index, m is the mass of the given fiavor quark, and 
are the Minkowski Dirac matrices. The Lagrangian in Eq. (3.1) is invariant under 
the SU{3) rotation [15] 

ijj{x) i/j'{x) = Aip{x), i){x) i>'{x) = i>{x)K\ (3.2) 

where A e SU (3) is independent of x. The above transformation is called a global 
gauge transformation [15]. If A depends on x, that is, A(a;) is a function of x, the 
Lagrangian in Eq. (3.1) is not invariant under the gauge transformation 

i^ix) i/j'ix) = A(x)'ijj{x), i){x) i)'{x) = i){x)K{x)\ (3.3) 

where A(a;) is an element of 5't/(3). The above transformation is known as local 
gauge transformation [15]. We will modify the Lagrangian in Eq. (3.1), so it it 
invariant under a local gauge transformation. 
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For each space-time point x, there is a space allowed values of iIj{x) at that 
point. How do we compare 'il){x) on one vector space \4 with tl^iy) at an other vector 
space Vy. In the gauge theory of QCD it is the parallel transporter that connects 
a point on one vector space with the corresponding point on another vector space. 
The parallel transport is generated infinitesimally [17] along a smooth curve. First 
let Cyx be a smooth curve in space-time from point x to point y. Then we define a 
mapping from vector space to vector space V^, that is, 

U{Cyx) : K ^ Vy, (3.4) 

for U{Cyx) € SU{3) such that U{Cyx)ip{x) is an element of the vector space Vy. 
That is, vector U{Cyx)ip{x) is regarded as the vector ip{x) at the point x, parallel 
transported along the smooth curve Cyx to the point y [15]. We say UiCyx) is a 
"parallel transporter" [15] along the smooth curve Cyx- 

We can prove that [15, 17], under the local gauge transformation 

V'(x) ^ V'(^) = A(x)V'(x), i^{x) i)\x) = i){x)K\x), (3.5) 

^{y) ^ ^'{y) = Ky)^{y). Hy) ^ i^'iy) = i'{y)^\y), (3.6) 

a parallel transporter transforms as [15] 

U{Cyx) ^ U\Cyx) = K\y)U{Cyx)K{x). (3.7) 

Now we consider the straight smooth curve Cx+dx,x from point x to point x-\-dx. 

A" 

The Lie algebra of SU{3) is generated by -— , the eight Gell-Mann matrices satisfy 
the commutation relations 

[A«,A^] = 2i/„,,A^ (3.8) 

where /"''^ are the structure constants of SU (3). Hence, the corresponding infinites- 
imal parallel transporter can be written as [15] 

U{Cx+dx,x) = 1 + igAl{x)^dx'', (3.9) 

where 1 is the 3x3 unit matrix, g is the color coupling constant, which is just 
introduced conventionally, A"^{x) is the gluon gauge field, and a is the gluon color 
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index (a = 1, .... 8). The parallel transporter in Eq. (3.9) has a vector index fi that 
gets contracted between ^^(x) and dx'^. For notational simplicity, we define the 
gauge field An{x) in the Lie algebra of SU{S) [15] 

A,{x)^zgAl{x)^. (3.10) 

Then we can rewrite Eq. (3.9) as 

U{C,+dx,x)^l + A^{x)dx''. (3.11) 

Now we perform a longer distance parallel transport along smooth path Cy^, where 
the points x, Xi,X2, x-^ - ■ ■ , Xn, y form a series of infinitesimally separated points 
along the path Cy^x- Then the parallel transport of '0(x) to the end of the path (i.e., 
point y) in space-time is constructed by taking a series of infinitesimal transports 
along path Cyx [15, 18, 19]. That is, 

U{Cy,x)^{x) = [1 + A^{x){x>t-x^)] X [1 + A^{x,){x^2-x^^)] x • • • x 

[l + A,{xn-i){x^r.-<-i)\ X [l + A,{xn){y^-x^M{y), (3.12) 

here U{Cy^x) is the parallel transporter along smooth curve Cy^x- This product of a 
series of infinitesimal parallel transports, with x (the path's beginning) on the left 
and y (the path's end) on the right, results in Dyson's formula [20, 15] 

U{Cyx) = P exp I / A^dx^ \ , (3. 13) 

where the symbol P denotes the path ordering of the integrand. The parallel 
transporter U{Cyx) is called the Wilson line from point y to point x along path Cyx, 
and when the path Cyx is a closed loop, the trace of the parallel transporter U {Cyx) 
is referred to as a Wilson loop by convention. 

Since now the parallel transporter can relate different vector spaces (14), we can 
define the covariant derivative of ip{x) by the formula 

D^ij{x)dx^' = U{Cx+dx,x)ip{x + dx) - ij{x). (3.14) 

After some algebra, we arrive at [15] 

D^^d^ + A^. (3.15) 
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We can prove that, under the local gauge transformation, the covariant derivative 
transforms as [15] 

D^iPix) ^ D'^iP'ix) = A''{x)D^iP{x). (3.16) 

Next we introduce the field tensor F^,^ through the commutator of two covariant 
derivatives [15]: 

F^,{x) = [D^, D,] = d'^A^ix) - d^'A^ix) + [A^{x),A,{x)]. (3.17) 

Since the field tensor F^j^ is also an element in the Lie algebra of color SU{Z), it 
can be written in the form [15] 

F,,{x)=igF;M^, (3.18) 

where = di^Al — d'^A^ — gf°-^'^A^^Al is the gluon field strength tensor. We 
can prove that, under the local gauge transformation, the field strength tensor 
transforms as [15] 

F^,{x) ^ F;,(x) = K\x)F^,{x)K{x). (3.19) 

If we replace the ordinary four-derivative by the covariant derivative D^^, the 
resulting new Lagrangian [13] 

3 

= E4(^)(^7^^M - mf)il^,{x) (3.20) 

c=l 

is invariant under the local gauge transformations. 

Since we have already introduced the gauge field A^, we should add its own free 
Lagrangian to Eq. (3.20). This term is given by the Yang- Mills Lagrangian (or pure 

gauge Lagrangian) [13] 

^^yM = -\FrF;,, (3.21) 

which is obviously invariant under the local gauge transformations. Hence, we 
finally arrive at the full gauge invariant classical Minkowski Lagrangian for QCD 
(here we consider the case of various flavors) 

1 3 ^ 

C = --Fr^;, + E E4^(^^ - (3.22) 

^ /=1 c=l 

where Ip — 7j^I?^, / is the quark flavor index, m/ is the mass of the quark with 
flavor /, ipl and i/S^ are the spin quark fields, is the number of the fermion 
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flavors, and 7^^ are the Minkowski Dirac matrices. For brevity, we have omitted 
the quark spin index. 

3.1.2 Continuum QCD formalism 

We calculate the physical vacuum expectation values of any observables F{U, t/j, -0) 
which are function of the field variables. For QCD in the Minkowski space [16], the 
expectation values can be written as a path integral ^ 







F{U, V', -0) I 0) = ^ J VUV^IjVxIj F{U, ip, -0) e'^^^ , (3.23) 



where -^and-^ are Grassmann variables ^, |0) stands for the ground state of the 
system (or physical vacuum states), and 

Z ^ jvUVi^Vi) e'^^. (3.24) 
Here Sj^4 — j d^x L is the action of QCD, L is described by Eq. (3.22), and 

VIJ = n^^M(^) (3-25) 

where V'andVJ are Grassmann variables, and a,^ are the indices including color, 
flavor, and spin. 

Hence, in Minkowski space the paths are weighted by a function e*'^ (5* is 
action) [16], which is oscillating. Therefore, it is not suitable for numerical calcu- 
lations. For a numerical simulation, it is very convenient to calculate in imaginary 
time (i.e., t — > — ir), where the variable r is called "Euchdean time." This 
means that we change the Minkowski spacetime metric to Euclidean spacetime 



■"^The path integral (PI) was introduced by Feynman in quantum mechanics in Ref. [21], For 
example, we replace the matrix element (a; |e~*^*| y) by an infinite-dimensional integration over 
all classical paths, which were weighted by the exponential of i times the classical action J Dxe^^, 
where the \y) is the eigenstate of Hamiltonian operator H, and S is the classical action. Hence, 
the Hamiltonian operator has been removed [15, 16]. 

^In Ref. [15] the authors give the detailed discussion of the properties of Grassmann variables. 
Please see Ref. [15] for details. 
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metric, then use the Minkowski Lagrangian in Eq. (3.22) for QCD to create the 
corresponding classical Euclidean action for QCD. First of all, we should understand 
the difference between these two metrics. We know that the Minkowski spacetime 
is a four-dimensional space (i.e., Xo,Xi,X2,X3), whose metric is g^^. ^ Hence, the 
Minkowski scalar product is 

x*y = g^^x'^y" = - - xV - ^V- (3-27) 

Euchdean spacetime (or Cartesian space) is the space of 4-tuples of real numbers 
(i.e., xi,X2,X3,X4), where X4^ (or t) is time and xi,X2,xs is space. In Euclidean 
spacetime, we use the metric 5fj_i,. Hence, the Euclidean scalar product is 

x-y = S^.x^y" = xY + + + ^Y- (3-28) 

In fact, Minkowski spacetime can be considered to have a Euclidean metric but 
with imaginary time x° = —icx^, where c is the speed of light. By convention in 
high energy physics, we choose c=l (so, x° = —ix^). 

Now we first make following substitutions: xq —1x4^ (or t —ir) and Aq — )■ 
—iA^, where t = Xq, and X4 = t. Second we adopt the Euclidean Dirac matrices 
7''. They are related to the Minkowski Dirac matrices through [15] 

7^if = il^ /^ = 1,2,3 (3.29) 

1% ^ 1^ // = 0. (3.30) 

Therefore, we obtain the classical Euchdean action [16] for QCD ^ 

/■ f 1 3 _ \ 

S^jd'xl -F^^F;^ + g E ¥cMM \ , (3.31) 

where 

^The definition of is in Ref. [22]. 

^In principle wc can relate Minkowski space observable to Euclidean space through analytic 
continuation. This is a trivial for static quantities such as masses [15, 22, 23]. In Chapter 5, we 
will discuss some problems (e.g., the energy conservation) due to this continuation. 
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Mf = Ip + ruf (3.32) 
$ = I'^D,, (3.33) 

and Mj is the Dirac operator. Then Eq. (3.24) can be rewritten as [16] 

Z ^ j VUVipVip e-^. (3.34) 

Therefore, Eq. (3.23) can be rewritten as [16] 







1 /■ - - 

F(C/,V,'0)|o) = - j VUVi^Vi^ F{U,i^,i)) e-^, (3.35) 



where S is the action in Euchdean space. We note that the weighting factor in 
path integral now is . 

Now we associate the path integral formula of the partition function described 
by Eq. (3.34) and Eq. (3.35) with the ones of quantum statistical mechanics. In 
quantum statistical mechanics, we define the partition function as 

Z{t) = Tre-^-^. (3.36) 

In terms of the eigenvectors and eigenvalues of the operator H which satisfies 
H\n) — En\n), where H is the Hamiltonian, r is the "Euclidean time." If we 
consider the following equation [15] 



Tr(e-^^F) = e"^"^ (n |F| n) (3.37) 

n=0 

Z(r) = Tr(e-"^) = £ 6"^"^ (3.38) 

We obtain that the thermal expectation values of any observables (F) through [15] 

Tr { Fe-^") 

(F) = lim — ^ ' . (3.39) 
\ I Z{t) ^ ' 

Hence, many calculations in QCD take on forms familiar in classical statistical 
mechanics. In Euclidean space the definition of the partition function described by 
Eq. (3.34) and the definition of the expectation value described by Eq. (3.35) are 
consistent with what we define in quantum statistical mechanics. Therefore, from 
now on we can reasonably call the Z in Eq. (3.34) "partition function." 
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We usually separate the action S into two parts 

S = S, + Sf. (3.40) 

The first term, Sg, depends on only the gauge fields. We call this term the pure 
gauge action (or Yang-Mills action) 

Sg^J d^x-^F^^F,,. (3.41) 

The second term, Sf, depends on both the gluon and fermion fields. We refer to 
this term as fermion action, 

Sf = j d'^x'^Mi/j, (3.42) 

where, for brevity, the flavor index and color index have been suppressed. 

For Uf degenerate fermion flavors, the QCD partition function in Euclidean space 
time is described by [16] 

Z ^ j VUV'^V^e-^^'^^i\ (3.43) 
Integrating the fermion fleld for degenerate flavors gives [16] ^ 

j Vi/jVi^e-^f = {detM)''f. (3.45) 

Thus 

Z ^ J VU{deiMYf e-^^ ^ J VU e'^-ff , (3.46) 
where the effective gauge action S^ff is denoted as [16] [24] 

Seff = Sg-Uf In (det M) . (3.47) 

The quark determinant detM generates quark loops. In the quenched approxi- 
mation, we set this determinant to a constant. In effect, in this approximation 



^For the n/ nondegenerate fermion flavors, 



/-J 
VipVipe-^f = JJ detM/, (3.44) 



where / is the index of flavor, and Mf is Dirac operator for flavor /. 
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wc ignore the interaction between fermions and gauge fields. This approximation 
introduces a 5 — 15% error in the lattice calculations of the light hadron masses [25] . 
It greatly simplifies the computation. 

The expectation value of any physical observable O that depends on the quark 
and gluon fields is [16] 



The above integral is a Feynman functional integral over all configurations of the 
gauge field. Usually the fermion fields are integrated out analytically, leading to a 
numerical path integral over only the gauge fields. 



In QED, a perturbative expansion gives quite accurate results for many physical 
problems. At very high energies, a perturbative expansion in QCD also gives quite 
accurate results for some physical systems. However in a low energy physical sys- 
tem, perturbation theory in QCD is not suitable, due to a large coupling constant. 
To overcome this, numerical techniques (namely, Monte Carlo techniques and a 
lattice discretization) were proposed to study QCD nonperturbatively. 

Numerical techniques require discretizing the continuum QCD theory on a Eu- 
clidean space time grid (a is lattice spacing). At the same time we should retain 
as many properties of the continuum QCD action as possible, and we require that 
the continuum theory be recovered in the physical continuum limit (i.e., a — > 0). 
Such an lattice QCD action is then used to generate a biased sample of gauge field 
configurations. Our observables are calculated on each of these configurations, and 
our observables are averaged on these configurations to get their expectation values. 

In this discretization, we optionally choose periodic boundary conditions for the 
gauge fields in space. That means that the first and last points along any straight 
line in the lattice are considered nearest neighbors (e.g., 0(xo) = (I){xn), where (f) is 
any function which depend on lattice site x. However in the time direction, periodic 
boundary conditions for the gauge fields are not optional. 



(O) 




-s 



-5 



(3.48) 



3.2 QCD on a lattice 
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In QCD perturbation theory, we encounter various divergent integrals. We deal 
with them by modifying the original integral. A "regulator" parameter characterizes 
the modification. The physical limit involves removing the regulator. The lattice 
is a natural regulator of ultraviolet divergence. In lattice regularization, we assume 

that our physical system is not a continuum space, but a discrete lattice. The lattice 

vr 

spacing a serves to limit the highest possible momentum (namely, — ) of particles 

a 

defined on the lattice. Due to this cutoff, all of the volume integrals of interest are 
sums: 



Hence, we never integrate the momentum to infinity. This shows that, in order to 
get the theory in the continuum, after performing the finite sums, we should take 
the limit a — >■ 0. 

In this section, we first give an overview of the pure gauge action in lattice 
QCD. Then we give three methods to incorporate the fermion, namely, the naive 
fermion, the Wilson fermion, and the staggered fermion. Finally, we discuss the 
taste interpretation of the staggered fermion. 



In order to explain confinement, in 1974, Wilson [26] proposed a lattice method. 
We now called it lattice QCD (LQCD). From then on, LQCD become a fundamental 
tool for the studies of QCD, and it is still in intensive use. The starting point for 
LQCD is to define fermion fields only on the sites of a four-dimensional grid with 
four spacetime dimensions of finite extent, N^, Ny, N^, Nf, and lattice spacing a. 
Here A^^, Ny, are space dimensions, and is the time dimension. 

A crucial point in his paper [26] is the introduction of the gauge link matrices that 
connect nearest neighbor lattice sites. Here we follow the procedures in Sec. 3.1.1 
and Sec. 3.2.2 in Ref. [15]. But in some places we choose different convention. 

For the lattice, let x be any point on the lattice, and lattice spacing be a. Then 
X + ap, is the neighboring point in the direction of the lattice axis /i — 1,2,3,4. 




X 



(3.49) 



3.2.1 Pure gauge action on lattice 
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where fi is the unit vector in the /i direction. The parallel transporter that moves 
along the straight path line from point x to point x + ap, is [15] 



We call Un{x) the link matrix linking site x to its nearest neighbor in the /i direction, 
X + ap,. A link matrix is an element of the gauge group SU (3) between pairs of 
nearest neighbor lattice points. The gauge field configuration is just the set of all 
such link matrices between pairs of the nearest neighbor sites. 

An important property in lattice QCD is that the trace of any product of link 
matrices around any closed loop is gauge invariant. This property provides a way 
to translate some gauge invariant quantities onto the lattice. For example, the 
smallest closed loop on the lattice is the unit square (see Fig. 3.1), which contains 
four points, x, x+ap, x+aD, x+ap+aP. The corresponding product of four link 
matrices [15], 



Up = U,{x)U^{x + aP)Ul{x + aP)Ul{x) = U,.U^+au,^Ul^ai.,.Ul^,, (3-52) 




(3.50) 



Here we make the approximation 




(3.51) 



is called the plaquette variable in the ji — v plane at lattice point x. 



X + au 



X + ap + au 



X X + ap 

Figure 3.1. A plaquette on a hypercubic lattice. 
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The Wilson gauge action, expressed in terms of the plaquettes, is [15] 

= E /5 {l - ^ (Tr + Tr ^7pi) I = ^ /3 1 1 - ^Re Tr [/p| , (3.53) 

where the sum over P stands for the sum over all plaquettes, including each 
plaquette with only one orientation [15], that is, 

E = E E • (3-54) 

P X l<//<t/<4 

6 

In Wilson action /3 — where g is the bare couphng constant. The Taylor 

9^ 

expansion of the link matrix Ufj,{x) is [15] 

U^{x) = 1 + aA^{x) + ^M^f- (3-55) 
The Taylor expansion of the gauge fields is [15] 

{x + afi) = (x) + aA^A^ (x) , (3.56) 
where Al is the forward lattice derivative. That means 

AiA^ (x) = -4.(x + aA)-^.(x) _ ^3_^^^ 

a 

If we use the well-known Campbell-Baker-Hausdorf identity [15] 
and the formula [15] 



F,u^igF;,,- (3.59) 



Y 

F^, = Aj^A^ix) - AiA^ix) + [A^{x),A,{x)] (3.60) 

Tr(F^^F^,) = -^Fr^;,, (3.61) 

we find, after some algebra, that as a — > 0, the action in Eq. (3.53) can be changed 
into [15] 

SG^\T.^'FrF^^, + 0{a'). (3.62) 

Therefore, we indeed reproduce continuum QCD action as a ^ 0, and it is the 
correct QCD action for the gauge fields in the lattice. 
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If we add some terms around larger loops than plaquettes, as in the Symanzik 
improved gauge action [27] , the O {a^) error in the gauge action can be eliminated. 
The gauge action used in this dissertation was a Symanzik improved action [28], 



Sg — 



imp 



x;ii<u 



(1+0.4805 a,) J2 



1 



0.03325 a, ^ C^A , 



x;fj,<iy<a 



(3.63) 

where P is the standard plaquettc in the fi, u plane, R is the real part of the trace of 
the ordered product of SU{3) link matrices along 1x2 rectangles (see Eq. (3.64)), 
C is the real part of the trace of the ordered product of SU (3) link matrices around 
1x1x1 paths (see Eq. (3.65)), 



R 



pi/ 



-ReTr 
3 



-ReTr 



(3.64) 



(3.65) 



and here Pimp — 
factor 



10 



41og(^o) 
3.0684 



. The symbol uq is the tadpole improvement 



uo = ( -Re Tr{P) 



(3.66) 



In lattice perturbation theory, there exists ultraviolet (UV) divergence from tadpole- 
type graphs [28]. This contribution can be partially removed by absorbing them in 
lattice coupling constants. The common method is to replace lattice gauge matrix 
U{x) by U{x)/uo [28]. In this dissertation, we also use this method. 

3.2.2 Fermion action on lattice 

The naive approach to formulate the fermion action is simply to replace the 
symmetric derivative in the Dirac operator with a symmetric difference. Here we 
follow the procedure and convention as in Rcf. [16]. We can change the continuum 
action in Eq. (3.42) into the naive lattice fermion action [16], that is. 



f 1 ^ - 

Sn = a'^J^i'^Yl U^{x)i>{x+afi) - Ul{x-afi)i'{x-afi) 



(3.67) 
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where the symmetric difference is denoted as [16] 

1 



Dn^ — 

" 2a 



Un{x)'i{j{x+ajl) — Uj^{x—ajl)'ijj(x—ajl) . (3.68) 



This is called the naive fermion action. 

Unfortunately, in the continuum limit, the naive action brings out an unphysical 
doubling of the fermion modes in each lattice direction. This mean that for every 
physical mode at momentum A;^ = 0, there is a degenerate mode at momentum 

TT TT 

kn = —, where — is the cutoff momentum of the lattice. Hence, in four space- 

a a 

time dimensions, there arc 16 degenerate modes [15]. Several methods have been 
proposed to solve this "doubling problem." 

Wilson added a second order derivative term to the naive action to break this 
degeneracy. This term includes a coefficient factor r, usually chosen to be 1. For 
the physical k/j, — mode, this term does not change the continuum limit, but this 

TT 

term suppresses the remaining 15 k/^ — — doubler modes [15]. The Wilson fermion 
action is [16] 

Sw ^ Sn- ^X!"^ + + Ufj,(x - afx)'il;{x - ajl) - 2V'(x)|. 

(3.69) 

After some algebra, we get [16] 

[ 4 _ 1 ^ - 

Sw = a'^^<{m+-)il!{x)^{x) + — J2'4'i^)['i'' + 7t^Ui^{x)^{x + ap,) 

For brevity, we define the hopping parameter k 

K = . (3.71) 

2am + 8r ^ ^ 

Also we rescale the fermion fields with 

i/j i)y/2k. (3.72) 
Then the action can be rewritten as 

Sw = Y^'>\^xMxy'^y, (3.73) 
xy 
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where M^y are the matrix elements of the fermion matrix [16], that is, 

M:cy = ^xy - f^Y^l^r + 7^] U ^{x)5^^y_ai, + [t - 7 J Ul{x - afi)5^_a[i,^- (3-74) 

Prom above equation, we can note that the diagonal term in M^y is unit. 

For the naive fermion action, due to the doubling problem, a single quark in fact 
describes 16 identical replicas of the quark. Each of these replication is called a 
"taste" [24]. 

Kogut and Susskind [29] developed another approach to solve the doubhng 
problem called the staggered fermion Kogut Susskind action. In this approach, 
the 16-fold doubling problem of the naive fermion action given in Eq. (3.67) can be 
reduced to a factor of four if we keep only a single fermion field component on each 
lattice site. This procedure reduces the number of tastes per quark fermion from 
16 to 4. Wc use another trick (discussed later in Sec. 3.3.2) to reduce the remaining 
4 fermions tastes to 1 in the determinant. 

To keep just one single fermion field component on each lattice site, we do spin- 
diagonalization [15] . Staggered fermions x ^'^^ defined by the local transformation. 

^{x)^T,x{x) i^{x)^x{x)Tl r, = 7^i7f7f7f. (3.75) 

where is a 4 ® 4 unitary matrix, which diagonalizes all the 7 matrices in the 
action. In terms of the staggered fermions x the quark action can be written as ^ 



f 1 4 

«5s = ^ < amtXt{.x)xt{.x)+- ^ a^^ 

X,t IJ,= 1 



Xt (x) Uf,{x)xt (x+afi) -xt (x+afi) Ul {x)xt (x) 



= J2xt{^)M'[Ul,yXt{y), (3.76) 

x,y,t 

where the t is the taste index of a given flavor, and the matrix M* is given by [16] 



M*[U]x,y = amtS^y + l- Y^^x,, 



(3.77) 



^ ii{x)^x,y—ajj, ^jiix OjjX)5x^y-\-afx • 

The 7^ matrices have been replaced by the phase factors 

= ^ ^ = 1,2,3,4). (3.78) 

^This formula is similar to Eq. (6.26) [16], but we use different notation and add a taste index. 
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From Eq. (3.76), and Eq. (3.77), wc note that the different spin components of x 
are decoupled, because the phase factor ax,n depends only on the site index and 
direction index and does not have a spin index. Here we drop the spin index on 
X, keeping only the color and taste degrees of freedom at each site. So this trick 
reduces the original 16-fold degeneracy of naive fermions by a factor of four. 

The fermion action used in this dissertation is the "Asqtad" action [9]. The link 
paths of "Asqtad" action are illustrated in Fig. 3.2. The first term is the simple 
link. The second term is three link staple, and the third term is five link staple. The 
second term and third term are to reduce the flavor symmetry [9] . The fourth term 

TT 

(or seven link staple) removes tree level couplings at the momentum — [9] . The fifth 

a 

term (or Lepage term) is the correction to the form factor at small momentum [9] . 
The sixth term (or Naik term, the third nearest neighbor) is to improve the quark 
dispersion relation [9]. If we set the zero momentum coupling to one and all the 
other momentum couplings to zero, and we introduce the tadpole factors ( — )^~^. 

Uq 

5 1 1 1 

where L is the length of path, we obtain the coefficients -, — ^, — Mn"^, Wn^, 

^ ^6 F > 8 16 ° 64 ° ' 384 ° ' 

"16^°^' ~24^°^' I'espectively. 

The fermion action in Eq. (3.76) is first term in the "Asqtad" action. Here is 
the second term: 

Ss2 = Y.{ (^mtXt{x)xt{x) + a^c/x Xt{x)U^{x)U^{x+P)Ul{x+fi)xt{x+afi) 

X,t I /i=l u^fj, '- 



-Xt {x+ail)U^{x+fi) Ul {x+u) Ul {x)xt {x) 
The rest are constructed similarly. 



(3.79) 




Figure 3.2. The simple link, three link staple, five link staple, seven link staple, 
Lepage term, and Naik term of Asqtad action. 
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3.2.3 Tastes of Kogut-Susskind staggered fermions 
Due to the phase factors ax,fx the Kogut-Susskind staggered fermion action has 
translational invariance under shift by 2a. Therefore, in the continuum hmit, we 
can map a 2^ hypercube into a single lattice point [15]. For simplicity, we just 
discuss the free Kogut-Susskind Staggered fermions. 

In position space [15], we can easily understand the taste identification. For 
some purpose it is convenient to represent the 16 degrees of freedom in a hypercube 
by a new field q called a hypercubic field, which in the continuum limit represents 
the four spin components of four degenerate tastes. 

Let us first assume that the number of the points in a hypercubic lattice is even 
in every direction: — 2L'^ [ji — 1, 2, 3, 4). The coordinates of the lattice points 
lie in the interval < x^, < — 1. The coordinates of the hypercubic blocks [15] 
are denoted by 

= 2y^ + ^^, (3.80) 

where — 0, 1. Therefore, it is convenient to divide the lattice into the hypercubes 
identified by a 4-vector that is, y is defined on a lattice of spacing 2a. The set 
of 16 vectors {A} identify the points within the hypercube, and they are called 
hypercube vectors. The sum over lattice points can be written as 

E - EE- (3-81) 

X y A 

The new fields with definite taste are written as g^" and g^", that is, 

C = I ErA;aa X^'^V + ^), C = 7 Ex(22/ + ^) ^aa, (3-82) 

^ A ^ A 

where a is the Dirac spin index, a is the taste index, and 

Ta = 7^72^^73^^7^- (3.83) 

The inverse expressions are 



X{2y + ^) = T:( F^g, ), x(2?/ + A) = Tr( g,F^ ). 



(3.84) 
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Then the free staggered fermion action can be written in terms of the new fields as 



Ss = 4:mqY^qy{l ® l)qy + Aj^Qy 7^ ® - 75 ® {^^,^5fD^ 



Qy, (3.85) 



y y,iJ- 
where we define two different lattice derivatives on the block lattice, 

= '^{fy+p, ~ fy-p.) ^ ^^ixfy 

= \Uy+P + fy-p-Vy)^a^dlfy. (3.86) 

In the continuum limit, the first term tends to the usual kinetic term. The second 
term is a lattice artifact of order 0(a), which disappears in continuum limit [15]. 
At finite lattice spacing a, the 16 degrees of staggered fermions in the hypercube 
are a mixture of spin and taste. 

We can write the general bilinear operator as 

OsF = Uls <H) ^T)qy = E X(Z/ + A) IV(r^75rBer) xiv + B), (3.87) 

A,B 

where = ^t- The matrices 75 determine the angular momentum and parity of 
the operator [30] and the matrices label the taste. Each of these can be one of 
the standard 16 elements of the Dirac Clifford algebra. 

The notation 7s stands for the direct product of the matrix 75 and 
where the former acts on the Dirac indices, and the later acts on the taste indices. 
The bilinear operators are identified as follows [31] [30] 

P = q{.l^®iT)q 

S = q{l(^^T)q, (3.88) 

where is a vector, ^4^ is an axial vector, P is a pseudoscalar, and 5" is a scalar 
meson [32]. 

We use these operators as interpolating operators to create and annihilate mesons. 
In the continuum, the mesons are in the taste SU (4) representations 1 and 15, while 
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on the lattice the action has a reduced symmetry and these representations break 
into many smaller representations [30]. For example, the continuum 15-plet of 
pion mesons breaks into seven lattice representations (four 3-dimensional and three 
1-dimensional) on the lattice. 

15 (75 (g) 7i) e (75 (g) 7,74) e (75 ® 7^75) ffi (75 <S) 7^7475) ® 
(75 ® 74) ® (75 75) ® (75 7475), 

1 ^ (75 (3.89) 

In the first hne of Eq. (3.89), i — 1,2, 3, the representations are three-dimensional. 
The other representations are one-dimensional [30]. Hence, when we measure the 
mass of the pion meson on lattice, there exists many different types of the mass for 
pion meson. In Chapter 8, we can see this in detail. 

3.3 Simulation algorithm 

In this dissertation, we use the configurations which were generated by the 
hybrid- molecular dynamics "R algorithm" [11]. In order to explain this algorithm 
clearly, we first introduce the process of generating gauge configurations. Then we 
discuss the hybrid-molecular dynamics "R algorithm." 

3.3.1 Numerical simulation and upadating process 

One of the main aims of numerical simulations in lattice QCD (LQCD) is to 
calculate the expectation values of some observables F[0] which are the functions 
of the field variables [0] = {(pxa}, where (pxa is a real field component with index 
a at lattice site x, where the index a includes flavor, color, spin [15]. For brevity 
we often omit a. For simplicity we assume that the variables for the fermions are 
integrated out analytically. 

im) = IT—, (3.90) 

/ [d0] e-^l-^l 

where e"'^''^] is an action suitable for simulation. The above equation involves a 
multidimensional integral over all the field variables [(f)] on the lattice. But the 
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number of the integration variables in [d(p] = d4>xa is typically very large (in this 

dissertation, it is over 10^). Hence, a Monte Carlo integration is used to calculate 
the integral efficiently [15]. 

For a given lattice with field variables [(f)], a configuration (i.e., [(f)]) is just one 
point in the space (or domain) of the field variables. Fortunately, due to the factor 
g-5'[</>]^ only those paths for which the action S[(f)] has the vahics close to minimum 
dominate the path integrals [15]. Hence, the efficient Monte Carlo simulation uses 
an importance sampling [15] [16], which generates a series of configurations with the 
density W[(f>] proportional to the Boltzmann factor e"'^'''''] [15]. That is, 

W[(f)] oc e-^l"^]. (3.91) 

We can call this series of configurations an ensemble of configurations. Then we 
can use these configurations to calculate the integral like Eq. (3.90) [15]. 

The configurations are generated {[0„], 1 < n < N} in a sequential manner 
by creating from one configuration to the previous configuration, [(f>] — >■ [(f)'], with 
a given transition probability. The transition probability, P{[(f>] — [(/)']), specifies 
the probability that the configuration becomes to a configuration 0'. The above 
procedure of creating configurations is called "updating" [15]. For the numerical 
simulation, we have a variety of schemes to get from one configuration to others. 

As in statistical mechanics, we want the density W [0] of the configuration. We 
also set the equilibrium density, Wc [(f>], to be proportional to the e''^^'^\ where S[(l)] 
is the action for a configuration. To achieve this result, a single step causes the 
density W to evolve as follow [15] 

W [(!>'] = Y.Pm^[(f>])W[(f>]. (3.92) 
The transition probability P satisfies [15] 



m 



(3.93) 



P([0']^[0]) > 0. 
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The normalization condition for the density W is [15] 

J2W [<!)'] = 1, (3.95) 
W] 

and we require the detailed balance condition [15] 

p m ^ [cf>]) [0] = p m ^ m m . (3.96) 

Different choices of the transition probability -P([0] — > [(f)']) generate different 
ensembles. Here we discuss one common updating algorithm: the Metropolis 
method [15] [16]. Briefly, the Metropohs proceeds as follows [15] [16]: we flrst 
randomly choose a possible new configuration (i.e., [(f)']) from possible new 
configurations, then calculate the new action (5'[0']). If the new action is less 
than the original action ( 5" [0]), we accept the configuration [(f)']. If the new action 

S[(f>'] is greater than the original action S[(f>], we accept the configuration [(f>'] with 

e-si4''] 

probability — ^ttt. We can simplify this algorithm for the transition probability 

P{[(f)] [(f)']) as follows [15] [16] 

Pm ^ [(f)]) = AT-imin |l, ^1 . (3.97) 

This method results in Wc [</>] oc e"'^'''^' , which means that the probability that any 
particular gauge configuration is included in the sample (or an ensemble of gauge 
configuration) is proportional to e~^'^'^\ 

With an ensemble of the sample gauge configurations, the average value of any 
operator (O) is simply calculated according to its average [16] 

— 1 ^ 1 

^ = j;^E<^[^n] + 0{-^). (3.98) 

If we have a large number of the sample configurations (i.e., N large), the average 
value that we measure is very close to the operator expectation value (O) . 
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3.3.2 Hybrid-molecular dynamics "R algorithm" 
The Metropolis method is simple. However, it usually updates the gauge link 
matrix locally (i.e., each step updates one gauge link matrix). Hence, we need a 
more efficient algorithm to update the gauge link matrix globally (i.e., each step up- 
dates the entire lattice) . The hybrid-molecular dynamics (HMD) "R algorithm" [11] 
is a example. It was used to generate the configurations for this dissertation. 

The basic idea of the HMD "R algorithm" can be stated as follow: by introducing 
a set of conjugate momenta H for the gauge link matrices U, the Euclidean path 
integral associated with QCD theory can be rewritten in the form of a partition 
function for a classical statistical mechanical system in space-time, and the motion 
of gauge link matrix (U) is governed by a classical Hamiltonian in the simulation 
time. Therefore, the configurations are generated in a deterministic way. 

In Sec. 3.3.1, we assume that the variables for the fermions are integrated out 
analytically. In some cases, we cannot analytically integrate out the variables 
for fermions. Hence, we consider the general case. The starting point for this 
algorithm [11] is to consider the following path integral for the partition function 
for Nf degenerate staggered fermion fiavors in Eq. (3.46). ^ 

/Nf 
[det{M{U))]~ 

= JvUe-^^ff, (3.100) 

where, as we explain below, 

Seff = Sg-^]ndet[M{UyM{U)]. (3.101) 

^In Sec. 3.2.2, when we discuss Kogut-Susskind fermions, we keep only a single fermion field 
component on each lattice site. This procedure reduces the number of tastes from 16 to 4. The 

factor - (fourth root trick [33]) is used to reduce the remaining four fermions tastes to one taste 

in the determinant. In this dissertation, we use the configurations [10] with three fiavors: two 
light quarks {u, d), and a strange (s) quark with = rrid mg. Then Eq. (3.100) becomes 

Z = j VUe-^^ [det(M„(?7))]5 [det{Ms{U))]^ , (3.99) 

where Mu[U]x,y is the quark matrix for uoi d quark, and Ms\U]x,y is the quark matrix for s quark. 
Hence, only two masses of KS quarks are included (i.e., the square root oi u ov d determinant, 
and the fourth root of s determinant). 



37 



Here Sg{U) is the pure gauge action on lattice, and the quark matrix M[U]x,y is 
given by Eq. (3.77). From Eq. (3.77), we obtain det[M(C/)] = det[M(C/)t]. The 
factor of - comes from the fact that det[M(C/)] stands for four tastes per quark 
fermion. Here, in Eq. (3.101), we considered the fact that the determinants of the 
submatrices on the even sites oi M{UyM{U) and the submatrices on the odd sites 
of M{UyM{U) are equal [34]. 

In order for the matrix Ux,fj, still to be an element of SU{3), the molecular 
dynamics equation ^ of the motion for Ux^n at site x should obey [11] 

UxAt)-^Hx^^ii)UxAi), (3-102) 

where Hx^t) is a traceless Hermitian matrix, and Ux,n{t) is the derivative of Ux,ix{t) 
with respect to the simulation time {t). Hence, as we promised above, we introduce 
the auxiliary field H conjugate to the gauge field A into the partition function [11] 

= J VUVHe-'^, (3.103) 
leading to the effective Hamiltonian [11] 

-H = \t^{H'') + Sg{U) - ^Trln[M(C/)tM(C/)], (3.104) 



andTr(i/2) = ^Tr(i/,%). 

X,IJ, 

Now we briefly illustrate the algorithm to generate a set of gauge field configu- 
rations with a probability distribution proportional to e"^, that is, 

P[U,H](xe-^. (3.105) 

To this end, we adopt two types of updating steps. First, the H field is initiahzed 
by the heat bath method [35] [36], that is, we construct a random H field [11] 

Hx,, = Khl^. (3.106) 
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The general Hamiltonian equations of motion are 



TTi — 



where 0, is field variable, ttj is the momenta cannonically conjugate to (^ij, and H[(I),-k] is the 
hamiltonian. 
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by setting each h°;^ to a complex Gaussian random number with standard deviation 



< >= 1. Thus, the H field is first updated globally. 



Second, the U field is updated by the equation of the molecular dynamics. Prom 
Eq. (3.102), if we first update the Hx^n, then we can easily obtain U. Since the H^^n 
is updated globally, the U field is updated at the whole lattice. Because it is very 
difficult to calculate the inverse of the matrix M{UyM{U), it is standard to use 
a noisy estimator method [37] [36] to compute it. The integration of the equations 
of motion proceeds by the leapfrog method. It generates the sequence U{t + iAt) 
and H{t + iAt), where i — 1,2, ■ ■ ■ , N. we take about N — ^ molecular-dynamics 
steps. The result is called a molecular dynamics trajectory. If we repeat two types 
of updating steps (namely, the heat bath steps and molecular dynamics steps) many 
times, we obtain an ensemble of configurations. 

In hybrid-molecular dynamics "R algorithm," since P[U,H] oc e~^, the expec- 
tation value of an observable F[U] can be written in the form 

[[dU][dH] e-^ F[U] 

{F[U]) = (3.107) 

[dU][dH] e-^ 



Because F[U] does not depend on the momentum H, and j[dH] e ^ ca e ^"^^ . 
Eq. (3.107) can be rewritten in the form 

/ [dU] e-^^ff F[U] 

{F[U]) = J— , (3.108) 

/ [dU] e-^^ff 

which mean that the probability distribution P[U] oc e'^^ff . 

This algorithm is apphcable to any kind of fermions [11]. For example, in this 
dissertation, we use configurations [10] with three fiavors: two light quarks {u, d), 
and a strange (s) quark. In our simulation, the light quarks are degenerate in mass, 
and we choose the different mass for strange quark. Hence, we replace the effective 
Hamiltonian % described by Eq. (3.104) by 

n - ^Tr(i/2) + -]^i:r\n[Mi{U)^Mi{U)] - ^ Tr ln[M,([/)tM,([/)], (3.109) 

where Mi{U) is the quark matrix for the light quark, and MgiU) is the quark 
matrix for the strange quark. The calculation is similar to what we discussed 
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before. The only difference is that we generate two gaussian random vectors, and 
use two conjugate-gradient ^ calculations (i.e., one for the light quark, one for the 
strange quark) for each molecular dynamics step. 

Prom the above algorithm, we can see that there are two sources of systematic 
errors. The first error is due to the nonzero time step (i.e.. At 7^ 0) in the 
integration, which results in the quadratic dependence of the error [11] (i.e., the 
error is of the order (At)^). This is because that we integrate the molecular 
dynamics equation (i.e., Eq. (3.102), Eq. (46) and Eq. (47) in Ref. [11]) with an 
error of order (Ai)^ at each time step. The second error comes from the conjugate 
gradient residual [11] (i.e., the required accuracy of conjugate gradient calculation), 
which we choose in the conjugate gradient calculation. If we choose the conjugate 
gradient residual to be small, we will take too much CPU time. Otherwise, if we 
prefer it to be the large, the accuracy is not what we want. So in the simulation, 
we choose a suitable value for conjugate gradient residual. 

Another error comes from our measurements. Because the conjugate gradient 
calculation consumes too much CPU time, we measure the potential (or the spec- 
trum) at the physical time separation St. In our simulation, we make a measurement 
at intervals of six simulation time units (namely, 5t = 6 At). So our measurements 
are not entirely statistically independent (i.e., our measurement is autocorrelated). 
We can reduce the effect of this autocorrelation of measurements by blocking several 
consecutive measurements before fitting [38]. In the latter chapter, when we fit the 
mass of the /o meson, we will consider the effect of autocorrelation. 



^In mathematics, the conjugate gradient is an algorithm for the numerical solution of linear 
systems. For example, we want to solve the following linear system 

Ax = b, 

where A is the n x n matrix, which is symmetric and positive definite, 6 is 1 x n matrix, and x 
is the solution. The conjugate gradient method is an iterative method. In this dissertation, we 
solve many systems like this. 
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LEADING BORN-OPPENHEIMER (LBO) 

The heavy hybrid mesons can be investigated directly by the numerical simula- 
tion, and by the Born-Oppenheimer (BO) approximation, which was proposed for 
the study of the hybrids states in Refs. [7, 8]. Here we focus on the leading order 
in the BO approximation, and ignore the higher-order terms. We use the Leading 
Born-Oppenheimer (LBO) approximation to study the spectroscopy of the heavy 
hybrid mesons (QQg), which consist of a heavy quark (Q) and a heavy antiquark 
(Q) plus a gluon. We are particularly interested in the spin-exotic hybrid mesons 
with J^'-' quantum numbers, which is not permitted in the quark model. The 
lightest one is believed to be with J^'-^ — 1"+ [5], which does not mingle with a 
nonhybrid QQ bound state [4]. 

In this chapter, we use the LBO approximation to study the spectroscopy of the 
hybrid mesons. We first introduce some notations from molecular physics. Then we 
study the LBO approximation. Next we review the APE smearing algorithm, and 
introduce the static potentials. Finally, we calculate the hybrid quarkonia through 
the LBO approximation. 

4.1 Notation from molecular physics 

We borrow the standard notation from diatomic molecular physics to classify the 
static potential (i.e., gluonic energies). Here we describe the gluon wave function 
in analogy with the electron wave function [4, 5]. The symmetries of the gluons 
include: 

• : The magnitude (denoted by A) of eigenvalues of the projection of the 
projection Jg • n of the total angular momentum Jg of the gluon field onto the 
separation axis R of the quark Q and antiquark Q with unit vector n — H/R. 
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The states with A = 0,1,2,... are typically denoted by the capital Greek 
letters E, 11, A, . . ., respectively. 

• CP : The combined operations of the charge conjugation (C) and spatial 
inversion (P) about the midpoint between quark Q and antiquark Q is also 
a symmetry. Its eigenvalue is usually defined by rj. States which are even 
(77 = 1) or odd (77 = —1) under the above parity-charge-conjugation operation 
(CP) are denoted by the subscripts g, and u, respectively. 

• it! : For the E states (A = 0), there is an extra symmetry. The wave function 
of the E states which are even (odd) under reflection in any plane passing 
through the separation axis (i.e., the plane contains the separation axis) are 
specified by a superscript +(— )• 

For the states with A 7^ 0, the energy of the gluons is not changed by operation R. 
Hence, those states are degenerate (namely, A doubhng) . We can hst the low- lying 
levels of the gluon field as E+, E~, E+, E~, n^, n„, A^, A^, etc [5]. In this 
dissertation, we only consider the states E+ and n„. 

4.2 Born-Oppenheimer approximation 

The hydrogen molecule ion, H^, consists of a single electron in the Coulombic 
field of two protons (see Fig. 4.1). Here we make the approximation that the 
nuclei (i.e., two protons) are at rest. Since the nuclear motion is much slower than 
the motion of the electron, the electronic wave function depends on the nuclear 
positions. The Hamiltonian is 

where Ri and R2 are the distances to the electron from the respective protons, and 
e is the charge of the electron [39]. 

For each fixed R, we can use some methods (e.g., the trial wave function in 
the variational principle) [39] to solve for the motion of the electron, obtaining the 
ground state energy V{R) of the system of H} as a function of R. We solve the 
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Figure 4.1. The hydrogen molecule ion, H2 
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Schrodinger equation for the motion of the nuclei using V{R) as the interacting 
potential. Then we can recover information about the positions and motions (e.g., 
rotation, vibration) of the nulcei [39]. In molecular physics this method is called 
the BO approximation. 

The heavy hybrid mesons can be studied similarly with the BO approximation 
(see Fig. 4.2). The slow heavy quarks QQ correspond to the nuclei in the diatomic 
molecule while the fast gluons correspond to the electrons [7] . 

To use the BO method, we first consider the quark Q and antiquark Q as the 
spatially fixed color field sources and use lattice methods to determine the energy 
levels of the gluons as a function of the QQ separation R. Each of these energy 
levels defines a static potential Vqq{R) [5]. At the second step in BO method, 
the heavy quark motion is then reproduced by solving the Schrodinger equation in 
each of these static potentials [5]. Conventional quarkonia are based on the static 
potential (S+), and the hybrid quarkonium states should come from the excited 
potentials (i.e., n„, etc). 



gluon 



9 



Quark Antiquark 




Figure 4.2. The hybrid meson, QgQ 
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In this chapter, we present the resuhs for the spectrum [5] of the gluonic excita- 
tions in the presence of a static quark-antiquark pair. Using our potentials, we also 
determine the hybrid quarkonium spectrum. For comparison with experiments, 
we also present the conventional quarkonium spectrum, which comes from 
potential. Before we evaluate the static potentials, we fattened [9] the spatial 
gauge links by "APE smearing" [40] [41]. Hence, we first discuss the APE smearing. 

4.3 APE smearing 

To reduce the mixings of our operators with its excited states the gluonic lattice 
operators are constructed from the iteratively-smeared spatial links (i.e., the fat 
hnk). We used the APE smearing algorithm described in Refs. [40, 41] to fatten 
the spatial gauge hnk matrix L'^(x), // = 1,2,3. The APE smearing averages a 
link matrix U^{x) on the lattice with its nearest neighbors. This means that the 
smearing procedure is to replace the ordinary link matrix U^{x) with the weighted 
sum of the link and a times its four staples (i.e., three link paths), that is, 

U^.{x) ^ U'^ix) = {l-a)U,{x) + '^j2[u.i^)U^i^ + i)a)Ul{x + i:ia) 

+Ul{x - va)U^{x - i>a)C/^(x + /ia)|, (4.2) 

where a is the smearing fraction. In our dissertation, we choose smearing fraction 
g 

a — — in order to get a satisfactory ground state enhancement [42, 43, 44, 45]. 

We construct the new gauge hnk U'^{x) as follows: Step 1: Calculate the new 
gauge hnk U'^{x) by Eq. (4.2) and project back to SU{?>). Step 2: After we do the 
step 1 for every link, we replace the old gauge link C/^(a;) by the new gauge link 
U'^{x). Step 1 and Step 2 define a single smearing sweep. One such a smearing 
iteration step is shown in Fig. 4.3. By convention, the new gauge link U'^{x) is 
called a "fat link" [46]. In this dissertation, before we evaluate the correlators, we 
take 20 smearing sweeps. 
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Figure 4.3. Visualization of a smearing iteration. 

4.4 Computation of the static potentials 

The first step in the LBO approximation is to determine the energy levels of 
the gluons as a function of where heavy quark {Q) and antiquark {Q) are 
fixed at lattice sites with separation distance R apart. At this point in the LBO 
approximation, the quark Q and antiquark Q which provide the color field are 
considered to be approximately static [5]. 

It is possible to describe the spatial distribution of the color field between the 
static quarks QQ. There are many papers [3] that study the physical properties 
of the color fiux between static quarks. We are very interested in the transverse 
extent of the color fiux tube (see Fig. 4.4) and the physical properties of the color 
fields (i.e., electric field or magnetic field). The ground state energy of the color 
electric field between static quark sources can be interpreted as the origin of the 
confining intcrquark potential of the heavy QQ pairs. The hybrid QQ potentials 
which bind quark-antiquark {QQ) pairs into the heavy hybrid QQg molecules are 
defined by the excitation content of the confining color fiux tube [5] . 

In this section, we first explain Monte Carlo estimates of the correlators, which 




Figure 4.4. The distribution of color electric field between static color charges. 
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arc used to calculate the static potential of the and Hu respectively. Then we 
discuss the procedure to fit these correlators. The fitting method we use is the 
standard Levenberg-Marquardt [47] fitting procedure. 

4.4.1 Caa Correlator 

To determine the potential energy of the QQ system we first need the correlator 
from the (A) state to the E+ (A) state. On a starting time slice t, the quark 
Q and antiquark Q are fixed at lattice sites with separation distance of R whose 
midpoint is at space point r (see Fig. 4.5). The gluonic operator at time slice t is 
the ground state, which is the parallel transporter for a straight link path along 
the quark and antiquark Before we evaluate the Caa correlator, we use the 

APE smearing method to "fatten" spatial finks, which means, in effect, that the 
parallel transporter here is not the conventional parallel transporter using the link 
matrix, and it is replaced by an weighted average over the link paths connecting 
the points [9] (see details in Ref. [46]). The gluonic operator at the source is 

Osourceir,R,t) = S+(r,R,t) = U{r-j+t4,r+^+t4), (4.3) 

where 4 is the unit vector in the time direction. 

At the final time slice tf (i.e., tf = t + T), the gluonic operator at time slice 
tf is the parallel transporter, which is a straight link path along the quark Q and 




Figure 4.5. An illustration of the correlator from A to A. 
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antiquark Q. The gluonic lattice operator at the sink can be written similarly as 



C»5mfe(r,R,t+r) = E+(r,R,i+r) = U{v-^+tfl,v+^+tf4.). (4.4) 

Therefore, the timeshce correlator Caa{T, R) from the state to the state 
can be written as 



Caa{T, R) = ^ E (Tr [cLik(r, R, t+T) Osourceir, R, t) 

T.t 



1 



n 



T,t 



R ^ R x+T^/ R R 
C/(r-_+i^4, r+-+i^4)tC/(r--+M, r+-+M) 



^(Tr[E+(r,R,t+T)tE+(r,R,t) 



(4.5) 



where Q = N^NyNzNt is the total lattice volume, Nx,Ny, are space dimensions, 
and Nt is time dimension. 

A few algebraic steps relates the correlator to the energy Vj.+ {R) of the lowest 
state excited by the operator. 



1 



Caa{T,R) = -^Tr(0 E+(r,t + T)S+^(r,t) 



n 
1 







^(0 E+(r,R,t + T) 



A) A 



S+^(r,R,t) 



-V^+{R)T 



e ^« 



E+(r,R,0) 



A) A 



E+^(r,R,0) 



where 1^4) is the eigenstate of the operator E+, and 



Za = {0 S+(r,R,0) 



A 



(4.6) 



(4.7) 



We have assumed that the contributions of higher excited states are negligible. This 
approximation is valid for large time T. 

Since the gauge configurations are gauge fixed to the temporal gauge (^4(3;) = 0) 
before we estimate the correlator Caa, the time-link link variables are unit matrices. 
Therefore, the product of parallel transporters in Eq. (4.5) is equal to the parallel 
transporter U{Cr^t) along the closed rectangular loop Cr^t of the side-lengths R 
and T (see Fig. 4.5), that is, Caa{T, R) = W{Cr,t): the standard Wilson loop. 
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4.4.2 Chh Correlator 

In this section, we calculate the correlator from the Hu state to the states (i.e., 
H state to H state). As before, on a starting time slice t, the quark Q and antiquark 
Q are fixed at lattice sites with separation distance of R whose midpoint is at space 
point r (see Fig. 4.6). The gluonic operator at time shce t is the n„ lattice operator, 
which traditionally is the parallel transporter along a U-shapcd path □ minus the 
parallel transporter along the same U-shaped path U (i.e., f/(Cn) — U{Cu)) in the 
opposite transverse direction (see Fig. 4.7). The symbol □ denotes the U-shaped 
path from the quark Q to the antiquark Q whose transverse extent can be in the 
positive X OT y direction, and U denotes the U-shaped path from the quark Q to 
the antiquark Q whose transverse extent can be in the negative x or y direction. 
The gluonic lattice operator at the source can be written as 




Figure 4.6. An illustration of the correlator from H to H. 
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X or y 



Q 



Q 



Q 



Q 



Figure 4.7. The creation operator for the n„ state. 

At the final time shce tf (i.e., tf = t + T), the gluonic operator at time shce tf 
can be written as 

Osinkir, R, t+T) = n„(r, R,t + T). (4.10) 
Hence, the timeshce correlator can be written as 

Chh{T, R) = ^ E (Tr \0lnki^, R, t+T) Osource{r, R, t) 



1 



r,t 



= - E (Tr [n„(r, R, t+T) n„(r, R, t))] 

T.t 



(4.11) 



Following steps similar to Eq. (4.6), we calculate the energy of the lowest n„ state 
Vn„(-R) through 



Chh{T,R) = ^ETr(0|n,(r,t + T)n,t(r,t) 







n 



Y,{0\U^ir,R,t + T)\H) (h Uj{r,R,t) 



r,t 



J2 e-^nJ«)T ^0 |n^(r, R, 0)1 H) (h IlJiv, R, 0) 



r,t 



- p-^n„(i?)r 

where \H) is the eigenstate of the operator n„, and 

Zh = {0\nu{T, R,0)\H) 



(4.12) 



(4.13) 



We have assumed that the contributions of higher excited states are negligible. This 
approximation is valid for large time T. 

Since the gauge configurations are gauge fixed to the temporal gauge (A^i^x) = 0) 
before we estimate the correlator Chh, the time- link link variables are unit matrices. 
Hence, in our simulation, we actually measure four bent Wilson loops (see Fig. 4.8). 
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/Q /Q Q Q /Q 

Figure 4.8. The bent loop for n„ states. 

4.4.3 The results for the static potentials 

In this dissertation, we want to reduce the mixings of our operators with their 
higher lying excited states. According we take 20 APE smearing sweeps. The 
correlators of CAA(t, R) and CHuit, R) are then fitted by using a single exponential 
term. Our best fit values are obtained by using the correlated method. For the 
correlator of CAAit^R), our fitting range of minimum distance tmin to tmax of the 
source-sink separation is listed in Table 4.1, where the a is the lattice spacing. 
And the average least per number of degrees of freedom for different R is 
approximately close to one. We hst some of the fitting information in Table 4.1. 

For the correlator of CHH(t,R), our fitting range of minimum distance t-^m to 
^max of the source-sink separation is listed in Table 4.2. And the average P^r 

Table 4.1. Fitting information for static potential V^+. The first column is the 
separation distance, the second column, static potential, the third column, the 
statistical error of the corresponding static potential, and the remaining columns 
are P^i' number of degrees of freedom, and the time range for the chosen fit. 



R/a 




Error for c/Tv 




Fitting range 


1 


0.5396 


7.999040e-04 


0.2/3 


8-12 


2 


0.8340 


5.111520e-03 


5.3/3 


7-11 


3 


0.9927 


7.192000e-03 


1.0/3 


6-10 


4 


1.0996 


1.826420e-02 


0.7/3 


6-10 


5 


1.2234 


1.291270e-02 


0.9/3 


5-9 


6 


1.3007 


2.307760e-02 


2.2/3 


5-9 


7 


1.3912 


4.568450C-02 


0.2/3 


5-9 


8 


1.5070 


7.961440e-02 


3.8/3 


5-9 


9 


1.6254 


3.601540e-02 


3.5/3 


4-8 


10 


1.7338 


5.708190e-02 


6.6/3 


4-8 


11 


1.8307 


1.003580e-01 


2.3/3 


4-8 


12 


1.8838 


2.552650e-02 


2.5/3 


3-7 
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Table 4.2. Fitting information for static potential Vn„. The first column is the 
separation distance, the second column, static potential, the third column, the 
statistical error of the corresponding static potential, and the remaining columns 
are per number of degrees of freedom, and the time range for the chosen fit. 



Rja 




Error for aVn 




Fitting range 


1 


1.6602 


6.8816e-03 


1.0/3 


3-7 


2 


1.6430 


6.1839C-03 


0.7/3 


3-7 


3 


1.6442 


6.0712e-03 


3.3/3 


3-7 


4 


1.6586 


6.5216e-03 


4.3/3 


3-7 


5 


1.6996 


7.6446e-03 


2.0/3 


3-7 


6 


1.7279 


9.1389e-03 


5.0/3 


3-7 


7 


1.7976 


1.1532e-02 


4.9/3 


3-7 


8 


1.8317 


1.5096e-02 


0.9/3 


3-7 


9 


1.9373 


2.1053C-02 


1.2/3 


3-7 


10 


1.9626 


3.0135C-02 


2.1/3 


3-7 


11 


2.0732 


4.1658e-02 


3.3/3 


3-7 


12 


2.1580 


6.1574e-02 


3.4/3 


3-7 



number of degrees of freedom for different R is approximately close to one. We hst 
some of the fitting information in Table 4.2. 

4.4.4 Setting the scale 

This scale parameter r\ depends only weakly on the valence quark mass and can 
be calculated accurately [48] [38]. Hence, we prefer to use it to rcscalc our data. 
The definition of this parameter and computation of this quantity are shown in 
Ref . [48] [38] . The scale parameter r\ is determined through 



R 



2 



dR 



= 1.0. 



R=rx 



From Ref. [10], we obtain the following value for r\. 



(4.14) 



ri = 0.317(7) fm or 
•f^ = 621(14) MeV. 



(4.15) 
(4.16) 



From Table II in Ref. [10], we get 



— = 2.632(13). 
a 



(4.17) 
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Therefore, we obtain [10] 

= 1.634(45)GeV. (4.18) 

Our fitting results for the static ground potential V-^+{R) and first gluonic excita- 
tions potential Vu^{R) are shown in Fig. 4.9. The results are expressed in terms of 
the scale parameter ri. 



4.4.5 Fitting procedure 

The static ground potential Vj.+ (i?) is shown as the low part (the octagon one) 
in Fig. 4.9, the solid curve for Vj^+{R) potential in Fig. 4.9 is a fit to the data by 
using a Coulombic plus linear form and lattice Coulombic potential [49], that is. 



+ aR, 



(4.19) 



8 — 



6 — 



> 4 



"1 — I — I — r 



"1 — I — I — r 



T — I — I — r 



T — I — I — r 




Q I I I I I I I I I I I I I I I I I I I I l_ 

2 4 6 8 

Figure 4.9. The static ground potential Vj^+ (R) and first gluonic excitations poten- 
tial Vyi^{R) in terms of the scale parameter ri with respect to the quark- ant iquark 
separation R. 
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is the lattice Coulombic 



where Vq, Cc, and a are three fitting parameters, and 
potential, which is calculated in free field theory for the Symanzik improved gauge 
action [50] [51] [49] as 



/TV 
-IT 



^ (Pk 



i=i ^ ^ i=\ ^ 



(4.20) 



This lattice correction term is used to correct the rotational symmetry of the 
lattice [49]. In practice we employ it at distances less than 3a. The Vq can be 
interpreted as a self-energy. Since only the tree-level coefficients are needed, we 
use Heller's data for the lattice correction term ^ 



R 



. In Table 4.3, we hst some of 



them. 

Our fitting gives an acceptable result with a of 46.9 for 39 degrees of freedom, 
with 

riVo = 2.2731(80) 
ec = -0.3927(27) 



0.5870(39), 



where we choose all the data within the fitting range of r from 



R„ 



0.0 to 



Rn 



ri 



ri 



7.89688. Since nVo = 2.2731(80), we obtain the self energy Vo, that is, 

Vo = 1.4116(45) GeV. (4.21) 



Table 4.3. for separate distance R. The first column is the separation distance. 



the second column is 



, and the remaining column is ^. 



R 


1 

R 


1 

R 


1.000000000 


1.038940053 


1.000000000 


2.000000000 


0.512582839 


0.500000000 


3.000000000 


0.332863472 


0.333333333 


4.000000000 


0.247627405 


0.250000000 


5.000000000 


0.197306796 


0.200000000 


6.000000000 


0.163888673 


0.166666667 


7.000000000 


0.140048006 


0.142857143 


8.000000000 


0.122177437 


0.125000000 
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The first gluonic excitations potential Vn„(-R) is illustrated as the upper part in 

Fig. 4.9. From Fig. 4.9 we note that at the short distance R the force between QQ 

is repulsive for the n„ potential, but the force between QQ is attractive for the 

potential. This can be explained by the color factor [13]. The potential 

(V2+) is the singlet static potential in which the quarks QQ state is in the color 

singlet representation. The S+ potential behaves like Eq. (4.19). On the other 

hand the !!„ potential is a hybrid static potential in which the QQ state is in the 

color octet representation. We know that the color factor fc — —-^ [13] for color 

— 4 

octet configuration of QQ, and fc — - [13] for color singlet configuration of QQ. 

3 

At short distance i2, the potential Vu^ should behave like — ^ [52]. That is, 

Vu. = ^, R^O, (4.22) 

where ag is strong coupling constant. Eq. (4.19) and Eq. (4.22) can help explain 
why at short distance R, the force between QQ is weakly repulsive for the Vn„ 
potential, but the force between QQ is attractive for the potential. 

The sohd curve for the Vu„{R) potential in Fig. 4.9 is a fit to the data by using 
the equation 

1 
R 



+ Co + V^o + biR + b2R^ (4.23) 



where e^, Cq, bo, bi, and 62 are five fitting parameters. It gives an acceptable fit 
with a of 44 with 37 degrees of freedom, with 

Con = 2.018 ±0.211 

borj = 5.313 ±0.866 

b^rl = -0.761 ±0.101 

b^rf = 0.561 ±0.028 

= 0.036 ±0.017. (4.24) 

For all R, in Ref . [5] , the authors confirmed that the first-excited potential of the 
color fiux tube is the Uu potential. Hence, we expect that the lowest lying hybrid 
mesons should come from the Vn„(-R) potential. 
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4.5 Hybrid quarkonium 

The second step in the LBO approximation is to reproduce the heavy quark mo- 
tion by solving the Schrodinger equation with the potentials of Vj.+ {R) or Vu„{R)- 
For notational simplicity, we adopt the spherical coordinates. The reduced mass /i 
of the two-body system is 

^ Mq + Mq 2 ' ^ ' 

where Mq is the mass of the heavy quark Q, Mq is the mass of the heavy antiquark 
Q, and we consider Mq — Mq The reduced radial Schrodinger equation of this 
system is given by [5, 53] 



dR^ 



+ 2fi[E- V,s(R)] u(R) = 0, (4.26) 

where Kff(-R) = Vqq{R) + ^ is the effective potential [5], Lqq is the orbital 
angular momentum of the quark Q and antiquark Q (discussed below in detail) [5] , 
and 

u(R) = p(R) R, (4.27) 

where p{R) is the radial wave function, and Vqq{R) is either Vj.+ {R) or Vn„(-R) [5]. 
The total angular momentum of the hybrid meson {QQg) is denoted by [5] 

J = L + S, (4.28) 

where S is the sum of the spins of the quark Q and antiquark Q, and L is the total 
orbital angular momentum of the hybrid meson. The orbital factor L is given by [5] 

L = Lqq + J., (4.29) 

where is the total angular momentum of the gluons. We choose the vector n to 
be unit vector along the axis of the quark Q and antiquark Q (or z axis). For the 
gluon state with A, we write the mean value of as [53] 



(4.30) 
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Since classically, L^^ = R x Pqq, it should be perpendicular to the n = H/R on 
average (i.e., (J^qq • = ). Hence, we deduce that [53] ((L — J^) • n) = 0, or, 

(L ■ n) = (J, ■ n) = A. (4.31) 

In the LBO approximation, the eigenvalue L{L+1) of and the eigenvalue S{S+1) 
of are also good quantum numbers. Hence, the total orbital angular momentum 
of the quark Q and antiquark Q is written as [5] 

(%) = (L^> + (V> - 2(l . J,) (4.32) 
= L(L + 1) + (jj) - 2A2. (4.33) 

For the potential of Vj.+ , there are no gluon excitations, so <^J^^ = 0. For the Vu„ 
potential, we have A = 1, so (J^ • n) = 1 [5]. Hence, Jg can choose any integer 
number which is larger than 1 (namely, Jg = 1, 2, 3, ■ ■ •). From Eq. (4.33), we note 
that if the system of the hybrid meson have more than one gluon, the energy level 
will be higher than that of the system with just one gluon [5]. However, here we 
focus on the decay of the hghtest exotic hybrid state. Therefore, we consider the 
system with just one gluon state [5] (namely, ^J^^ = 2). 

We know that A is the ^-component of the angular momentum of the gluon. 
By convention if A is positive, we call it a left-handed gluon state (|left)) which is 
not an eigenstate of a PC. If A is negative, we call it a right-handed gluon state 
(I right)) which is also not a PC eigenstate. However, if as we discuss in more detail 
in Appendix B, the linear combinations of left-handed and right-handed gluons 
states are |left) -I- e |right) [5], where e = ±1. are PC eigenstates. For E+, since 
A = 0, we have e — 1. Let 77 = ±1 be the PC quantum number of the gluons, 
which correspond to the subscripts g{u). In Appendix A, we prove that the parity 
and charge conjugation of the hybrid meson are given in the form [5] 

P = e (-1)^+' (4.34) 
C = er^ (-1)^+5. (4.35) 

In general the potentials evaluated from Wilson loops in lattice simulations are 
built up by the same diagrams as the self-energy for the Green's function [54] . The 
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self energy Vq renormalizes the quark mass so we absorb it into the definition of 
Mq. The difference between Vq and cq is physically significant, however we should 
remove the unwanted self-energy contributions. In practice we fit our data for the 

potential to the form in Eq. (4.19). The constant Vq is our approximate estimate 
of the self-energy contributions. 

For the n„ state, we have A = 1, from Eq. (4.31), we obtain that L > 1. The 
spectra of the heavy hybrid meson (bbg) built on the potentials of the 1^2+ ^"^^ 
are shown in Fig. 4.10. The energy levels for each state are listed in Table 4.4. 

From Ref. [55], we know that in the heavy quark limit, the Hamiltonian operators 
for E+, n„ in our LBO approximation is 
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10.8 
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BB threshold 
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9.8 






9.6 
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2g+(p) n,(p) n,(D) 



Figure 4.10. Spin-averaged bb spectrum in the LBO approximations. Solid lines 
specify the experimental measurements. Short dotted, short dashed lines indicate 
the S and P state masses obtained by solving the Schrodinger equation with Vj.+ 

by using rf ^ = 0.621 GeV and Mbi = 3.88 GeV for the b quark mass, respectively. 
Long dashed, dot-dashed lines indicate the hybrid quarkonium states obtained from 
the Vn„ {L = 1,2) potentials respectively. 
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Table 4.4. Quarkonium and hybrid quarkonium as computed. The first column 
gives the first six energy levels for E+(5'), the second column, the first six energy 
levels for E+(P), the third column, the first four energy levels for n„(P), the fourth 
column, the first four energy levels for Ilu{D). 



KiS) GeV 


E|(P) GeV 


n„(P) GeV 


UuiD) GeV 


9.4603 


9.87656 


10.8283 


10.96 


10.0183 


10.2821 


11.1025 


11.2332 


10.4076 


10.6194 


11.3694 


11.4972 


10.7349 


10.9186 


11.628 


11.7525 


11.0269 


11.1924 






11.2951 


11.4475 







H^^ = 2M«3 + ^ + y^+(P) (4.36) 

2 

Hu^ = 2M,o + ^ + yn„(P), (4.37) 

where M^q is called the rest mass, and M^i is referred as the kinetic mass [55]. The 
heavy quark rest mass M^o was determined by matching the known mass of the 
T{1S). 

Mrns) = 2M60 + Eq, (4.38) 

where the Eq is the energy of the ground state in the potential of E+. From our 
LBO approximation, we obtain Eo = 0.1711 GeV. We know My^is) = 9.4603 GeV. 
Hence, we estimate 

Mfto^ 4.64 GeV. (4.39) 

The heavy quark mass is tuned by using a directly measurable physical 
quantity. A conventional choice is to use the experimentally known T spectrum, 
namely, the splitting of 2S-1S (m£;(2S-lS)) and the splitting of IP- IS (mB(lP-lS)). 
Prom the LBO approximation, we obtain the calculated splitting of 2S-1S (mc(2S- 
IS)) and the splitting of IP-IS (mc(lP-lS)). By adjusting the kinetic mass of the 
b quark {Mb), we make both the difference of splitting of 2S-1S (namely, mc(2S- 
IS) - mB(2S-lS)) and the difference of splitting of IP- IS (namely, mc(lP-lS) - 
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mE(lP-lS)) as small as possible. Then we obtain the kinetic mass M^. In practice 
we minimize this formula 

f{Mh) = ^J[mc{2S-lS) - m^(2S-lS)]' + [mc(lP-lS) - m^(lP-lS)]' (4.40) 

to obtain the mass of the b quark (Mb). In our simulation we estimate 

MM = 3.88GeV. (4.41) 

Prom Fig. 4.10, we note that below the BB threshold, our LBO results are 
consistent with the spin-averaged experimental measurements, and above the BB 
threshold, our LBO results are not in good agreement with the spin-averaged 
experimental measurements. This disagreement probably comes from the fact that 
above the BB threshold, the T(45') and T(55'), etc, can decay to BB [12] and 
our model does not take those decays into accounts. However due to their exotic 
properties and heavy quark kinematics, our hybrid states decay only weakly to 
BB (discussed in Chapter 5). Therefore, our hybrid energy levels are presumably 
insensitive to this effect. This disagreement also comes from the fact that in our 
LBO approximation we ignore the high order terms. 

Prom Pig. 4.10, we note that the mass of the lowest-lying hybrid from the Vh„ 
is approximate 10.83 GeV. In fact we can obtain the precise value of lowest-lying 
hybrid niH from the LBO approximation based on the potential of Vn„, namely, 

mi/ = 10.8283 GeV. (4.42) 

Above 11 GeV, the LBO approximation based on Vn^ predicts several hybrid states 
with approximately uniform energy level separations about 200 — 300 MeV. This 
is a distinct feature of the hybrid quarkonium spectrum. 

The reduced radial wave functions for the IS and IP conventional quarkonium 
and that of the lowest-lying n„ hybrid quarkonium state are shown in Pig. 4.11. 
Por notational convenience, we call the reduced radial wave function for the 15" 
conventional quarkonium state Urj^{R), the reduced radial wave function for the IP 
conventional quarkoniun state u^^^{R), and the reduced radial wave function for the 
lowest- lying 11^ hybrid quarkonium state uh{R)- 
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Figure 4.11. Reduced radial wave functions versus quark-antiquark separation R. 
The solid and dotted curves specify the reduced radial wave functions for the IS" 
and IP states of S+, respectively, and the short dashed curve specifies the reduced 
radial wave function of Uu hybrid state, where all the wave function are normalized. 

From Fig. 4.11 we note another distinct feature of the heavy hybrid states. It 
is the shape of the reduced radial wave function. The hybrid reduced radial wave 
functions all vanish at the origin [R = 0). This can be explained by the repulsive 
Coulombic term in the potential Vqq{R) and the special centrifugal term Ves{R) in 
the LBO Schrodinger equation [7]. As depicted in Fig. 4.11, we also note that the 
hybrid Schrodinger state have a broader reduced radial wave function than those of 
the conventional quark-antiquark states. In Chapter 5, we use these LBO results 
to develop the decay model for the hybrid state. In Chapter 8, we use these results 
to calculate the decay rate. 




CHAPTER 5 
HYBRID MESON DECAY CHANNEL 

Most hadronic states are unstable under strong interactions. A comprehensive 
understanding of the hadronic decay from QCD involves developing the physical 
models to deal with the unstable hadrons and to extract their hadronic transition 
matrix elements. The physical models for studying the strong decays of the hadrons 
have been studied for a long time [56] [57]. Here we concentrate on the hadronic 
decay in the static limit, and employ lattice methods to calculate its transition 
matrix elements. 

For the decays of the hybrid mesons, we study the creation of a light quark- 
antiquark state from the gluonic field of the hybrid meson, then in the heavy quark 
limit we discuss its open decay channels with the emphasis on the XbS channel, 
where 5" is a scalar meson. This scalar channel (1""*") leads to tt+tt", 7r°7r°, or 
K+K-, K+K-. 

5.1 Hybrid states on the lattice 

In Chapter 4, we studied the properties of the physical hybrid states by solving 
the Schrodinger equation with these static potentials in the LBO approximation. 
Prom Fig. 4.10, we can see that the ground state of the color flux is in the state 
of E+. The first excited state is confirmed to be in the state of in Ref. [5], and 
we show this result in Fig. 4.10. We found that the mass of a hghtest hybrid spin 
exotic with J-^^ = 1~+ is at ruH — 10.83 GeV, about 1.37 GeV heavier than the 
mass of the T(IS'). Above this energy level there are many hybrid states, including 
the spin exotic states. 
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5.2 Hybrid meson decay model 

In this section, we give a brief review of the possible decay channels of a hybrid 
meson through the creation of a light quark- antiquark pair (qq) in the heavy quark 
limit. The jointed symmetry of the qq and the final gluonic field is identical with 
the symmetry of the initial representation of the gluonic field [4] . 

For a spin-exotic hybrid meson which contains an excited gluonic field n„ with 
A = 1 about the separation axis between QQ, we obtain, from Eq. (4.31), that 
the orbital angular momentum of the QQ pair should be greater than or equal to 
one (L > 1) [4]. To make the problem simple, we assume that the orbital angular 
momentum of the QQ is in the P-wave (namely, L — 1). 

In practice, the n„ state corresponds to an excited gluonic field with L-^^ = l+~ 
or L^'-'' = [58]. With the consideration of the spin of the heavy quark QQ, 
we obtain a set of eight degenerate states [4]. From the discussion in Chapter 2, 
we noted that, J-^*^ — 1"+, 0+" and 2+~ are the spin-exotic hybrid state [58]. 
Therefore, the spin-exotic hybrid meson has the heavy quark- antiquark pair in a 
spin triplet {S — 1). In any strong decay this is also conserved. 

In Chapter 4, we estimated the energy of the low hybrid state, namely. Eh — 
10.83 GeV. Therefore, if only considering the energy conservation. The spin-exotic 
hybrid meson with quantum number J^'" = have open decay channels: BB, 
BW, B*W, ribT}, rjbT]', XbV: XbS, T{ls)uj and T(ls)0 [4]. 

Page proposed one non-relativistic symmetrization selection rule [59], namely 
that the hybrid decay does not lead to two identical states: 

Hy^X + Y. (5.1) 

That is, the decay is suppressed if X and Y have an identical internal structure 
in all respects, as well as an identical radial and gluonic excitation (namely, each 
meson has orbital angular momentum L — (S-wave)). Prom this rule we can first 
exclude the decays to BB, BB* and B*B*. 

For the hybrid meson with J^'^ = 1"+, we consider its decay to a {Qq){Qq)- If 
each heavy-light meson Qq has L = 0, it is impossible, because the spin triplet 



62 



state has even CP, and the spin singlet state has A = [4]. If one or two 
heavy-hght mesons have a nonzero orbital excitation(L > 1), it is permitted from 
the consideration of the symmetry, but not permitted energetically. Prom Ref. [12], 
we note that experimental observation mass of the excited (L=l) B meson is 
5.732 GeV, and the mass of nonexcited (L=0) B meson is 5.279 GeV. Therefore, 
the hybrid exotic meson with mass Eh = 10.83 cannot decay to these two B mesons 
if one or two B mesons have a nonzero orbital excitation (L > 1). 

On the other hand, the decays of a hybrid meson with J^*" — 1""*" to {QQ){qq) 
are allowed energetically [4], because the excitation energy in the gluonic field is 
enough to produce a light quark-antiquark qq pair in a fiavor singlet state [4], for 
example, r), ou and scalar (S) channels. However, as we have discussed, its orbital 
angular momentum is conserved. Hence, the decays to the quarkonia of QQ are 
only into a Xb state (P-wave state) in the static limit. 

For the decays from a hybrid meson with the H„ representation to a gluonic field 
with the representation, the qq should have a H„ representation with A = 1 and 
CP — —1. Keep in mind that this orbital wave function of the qq of separation 
axis between QQ has no spin component. Therefore, if the qq is the scalar meson 
with J^^ = 0++, its orbital wave function should have odd CP [4]. 



CHAPTER 6 



OUR CALCULATIONAL METHOD 

In Chapter 5 we discussed the possible decay channels of a hybrid meson in the 
heavy quark limit. The goal of this dissertation is to study these hadronic decays. 
In this dissertation, we consider the decay, 

H — >A + B, (6.1) 

where H is unstable spin-exotic hybrid state of the static heavy quark. The hybrid 
meson H is in the n„ representation. The A and B stand for two decay products. 
The A is a conventional quarkonium state Xh in the representation. The S is a 
0++7r7r state [4]. 

In Chapter 4, we explained in detail the correlators Caa and Chhi which we used 
for calculating the static potential Vj.+ {R), l^n„(-R) respectively. In this chapter, 
we discuss in detail the correlators Cbb and Chab, which are used to calculate the 
transition amplitude x and the decay rate T in Chapters 7 and 8. 

6.1 Cbb correlator 

In this section we calculate the correlator from B state (0"'""'" tttt state) to B 
state. We first give the definition of the Cbb correlator. Then we introduce two 
methods to calculate it. Because we use replica trick in our calculation for the Cbb 
correlator, we first briefly introduce the replica trick. 

6.1.1 Replica trick in lattice 

As we discussed in Chapter 3, for the KS staggered fermion, we reduce the 
number of the tastes per flavor from 16 to 4. To further reduce from 4 to 1, we 
used the fourth root trick. At the quark level, the fourth root trick can be achieved 
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by multiplying each dynamical sea quark loop by a factor of ^ [60]. In Ref. [60], 
the authors apply this method to staggered chiral perturbation theory (SxPT). 
However in practice, when we use this method to work at all orders in SxPT, we 
cannot make sure that every time we do it correctly at higher order in S^PT. A 
solution for this is to use the replica trick. ^ Wc can state this trick as [60] follows: 

1) Generalize the number of each quark flavor to contain rii quarks {i—u, d, s) (i.e., 
Uu for up quark, for down quark, Ug for strange quark). Therefore, we have 
n — rii total quark for all flavors u, d, s. 

i=u,d,s 

2) Calculate physical quantities in which we are interested (e.g., correlation func- 
tion) as analytic functions of the n^. 

3) In the end, set each rii — -, i — u,d,s. 

This trick automatically performs the transition from four tastes to one taste per 
flavor for staggered fermion at all orders. We use this method extensively. 

6.1.2 The operator for the Cbb correlator 

On a starting time slice t, the quark Q and antiquark Q are flxed at lattice sites 
with separation distance of R, whose midpoint is at space point r. Let s be the 

spatial point of the tttt states at the source time, and s' corresponding position at 
the sink time (sec Fig. 6.1). In order to make our tttt state to be O"*""*" flavor singlet, 
the fermionic lattice operator at the source can be written as 

where g is the index of the taste replica, is the number of the taste replicas, and 
a is the color index. At the final time slice tf {tf = t + T), the fermionic lattice 
operator at the sink can be written similarly as 

, . , , _^ 4(s^t/)^^^,(s^t/)+4(s^t/)4(s^^/) 

b,g' V r 

where g' is the index of the taste replica and b is the color index. 



^The replica trick is originally proposed in Refs. [61] to study the spin glass. Now it has been 
widely used in the theory of the disordered systems such as QCD [62]. 
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Figure 6.1. An illustration of the correlator from B to B. 

Therefore, the time shce correlator Cbb from B (namely, O"*""*" tttt) to B for a 
chosen momentum k can be evaluated by the formula 

CBB{T,k)^^ Y: {Tr[Ol^,{s',t + T) O s our ce{s, t)) e^"^'^^-^'^ 



s,s ,t 

1 1 

Q n 



E E E E t+TKis', t+T) K»(s, t)ui{s, t)) + 



t s',s a,b g^g' 

4(s', i+r)4(s', t+T) d»(s, i)ci«(s, t)) (6.4) 

where Vt is the total lattice volume. Here we consider that the u^d quarks are 
degenerate in mass. We can simplify the correlator of in Eq. (6.4) as 

Cbb{T, k) = i- ^ ^ E E|(4(«'' ^+T) u'^is, t)<(s, t)) + 

""-^ t s',s a,b g,g' I 

4(s',t+T)<,(s',t+T) <(s,t)d^(s,t)) |e*''-(^-^'). 



(6.5) 

After we perform the Wick contractions of the fermion fields, and sum over the 
index of the taste replica [15], we obtain 

Cbb{T, k) = ^ E E ( E {m^i>\^': t+T- s', t+T)M-^\s, t; s, t)) e^^<^-^'^ ] (6.6) 

" t s',s I a,b ) 
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2n, 



E E I E {Mba (s', t+T; s, t)M„V (s, t; s', t+T)) e^«^-(-^') | , 

t s',s V a,6 ) 



where M is the fermion matrix. Therefore, the correlator Cbb(T', k) contains two 
parts. The first term in Eq. (6.6) is the disconnected contribution part (i.e., DC 
part). We label this term C^§. 

CEEiT, k) = ^ E E { E {M-a\^^ t; s, t)M^,\s', t+T; s', t+T)) e^^ (^M. (6.7) 

" t s',s I a,b ) 

The second term in Eq. (6.6) is the connected contribution for correlator Cbb- We 
label this term C^b- 

CEE{T,\^) = ^Y.T.\T.{Mab\^^t;s',t+T)M^^\s\t+T;s,t))^'^^^^^^^^ (6.8) 

" t s',s I a,h ) 

6.1.3 Direct method 

In our lattice simulation, we use the direct method to calculate the connected 
part Cf^(r,k). Eq. (6.7) can be changed into 

CiBiT, k) = ^ E E |Trc [m-\s, t; s', t+T)M-\s', t+T; s, t)] e^''-^-^') |, (6.9) 

where Tr^ is trace over the color index. We know that 

M-\s, t; s', t+T) = [M-\s', t+T; s, ^)]^(-l)^'-^ (6.10) 

where s = {s,t),s' = (s',i+T), M~^{s;s') is 3 x 3 color matrix, and (—1)*'"^ = 
(^_iy'i+s'2+s'3-si-s2-s3+T _ Finally, we can rewrite our connected part as 

CEEiT,k) = ^j:j:li-l)^'-^e^^<^-^'^x 

Tr^[M-\s,t;s',t+T)M-\s,t;s',t+T)^] |. (6.11) 

We use the conjugate gradient method (CG ) to obtain the required matrix element 
of the inverse fermion matrix M(s, t; s', t+T) by solving the hnear system 

^ M„,6(s", t"; s', t+T)M,-^(s', t+T; s, t) = 5a,c5t",tSs",. (6.12) 

b 

for M^-^(s', t+T; s, t), where the a, 6, c are color indices. Since we obtain the matrix 
element of M~^{s' ,t+T;s,t), we can easily calculate C^^(r,k) by Eq. (6.11) . 
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6.1.4 Random source method 
With the random source we compute both the DC part and CC part for the 
correlator Cbb- This method involves relating the desired correlator to a correlator 
involving random color fields ^: 

t s',s I 

l^EEE [e(s,^)M-„}(s,i;r;„gC(r'o,g X 
1^ E E E [ef (s', t+T)M^,\s', t+T; r', t')C^(r', t') 



7]' b,b' r',t' 



EE 

t s',s 



gik-(s-s') ^ 



EE 

+/ _/ fi 



EE {M;^>is,t;v',X)M^,}{s',t+T-,v',t')) x 



a,a' b,b' 



^ E E (C* (s, ^)e(r;„ i'o)er(s', i+T)0'(r', t')) 



7V2 



S (6.13) 



where Nj. is the total number of random sources, ^ is the complex vector of the 
gaussian random numbers, a, a', b, h' are color indices, and v( are the indices of 
the random sources. In the MILC code, we simulate the ^ to satisfy 

(ReC.) = i- ^ ReC: = (Im^,) = i- ^ Im^: = 



r r 

^Wlf = \ (6.14) 



where Re^ is the real part of ^, and Im^ is the imagery part of ^. We can prove 

1 



AT, 



(6.15) 



■Trl2Q'i^)^biy)C'W)ib'iy') = 5xy5ab5a'b'5x'y' +5ab'5xy'5ba'5yx', (6.16) 

where a is color index, and 77 is the index of random source. 

For 77 ry', the correlator in Eq. (6.13) gives the disconnected part, that is. 



CBr(r,k),^,, = C^f (T,k) X (1 - ^) . 



(6.17) 
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For T] = T]', it gives both the disconnected part and connected part, that is, 

CE^iT, k),=,, = i- {cEEiT, k) + Cfg(T, k)) . (6.18) 

Therefore, we obtain both the disconnected part and the connected part of the cor- 
relator Cbb from this method. In practice, we compare the result of the connected 
part from this method with the result of the connected part from the direct method 
to justify our results from this method. Prom our calculation we find that if we 
choose the Nr large enough, the numbers from two methods agree. 
Let 

^'{^,t) ^ E E er(s.^)K-a'(s,i;r'o,i(,)C(r'o,i[)) (6.19) 

a,a' ro'jty 

a^'{s',t+T) ^ EEer(s',t+T)M,-,,^(s',t+T;r',t')e.^'(r',0- (6-20) 

6,6' r',t' 

Then Eq. (6.13) can be rewritten as 

Cfrmk) = ^^EEE(^ns,^)^^'(s',t+T)e^>^-(-')|. (6.21) 

ri,v' t s',s I J 

If we define the Fourier transform: a''{k,t) = ^a''{s,t)e~''^'^, then Eq. (6.13) can 

s 

be rewritten as 

C^r{T,\^) = ^4EE^''(-k,^)<T'''(k, W). (6.22) 

r t ri,ri' 

From the above discussion, we know that, for r} ^ r]', the above correlator gives the 
disconnected part and for rj — rj', it gives both disconnected and connected parts. 

6.2 Chab correlator 

In this section we show how to calculate the correlator from the 11^ state or H 
to our two final states A and B state. On a starting time slice the quark Q and 
antiquark Q are fixed on lattice sites with a separation distance R whose midpoint 
is at space point r (see Fig. 6.2). The gluonic operator at time slice t is the 
lattice operator as we described in Sec. 4.4.2. 

At the final time slice i+T, the gluonic lattice operator can be written as 

Osink{r, R, t+T) = E S+(r, R, t+T)w{s - r) x 
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Figure 6.2. An illustration of the correlator from H to AB. 



E 



u'^is, i+r)w"(s, i+r)+d"(s, i+r)(i"(s, t+T) 
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= ^/2 ^ S+(r, R, t+T)w{s-r)u''{s, t+T)u''{s, t+T), (6.23) 

s,a 

where we consider that the u, d quarks are degenerate in mass and their production 
amphtude are identical. The w(s— r) ^ is the wave function of ^"(s, t+T)u°'{s, t+T). 
Here we should note that wave function ^(s — r) is just our choice for the final 
state. We assume it is a real function (i.e., w* — w). We should choose a suitable 
wave function w(s— r) to make sure that the initial and final states are in same 
symmetries. With these definitions the time slice correlator can be written as 

Chab{T, R) = ^ E (Tr [0kik(r, R, t+T)Osource{r, R, t) 



y5 



E E (E ^''(s, t+T)u'^(s, t+T)\w(s-v) X 

r,t s \ a I 



'Iv[E+(r,R,i+r)n„(r,R,i) 
= ^i:i:Y.Maa\^^t+T-s,t+T)w{s-r) 



X 



r,t s a 



^In Chapter 5, we discussed the CP symmetry of the wave function w{s — r). 



Tr 



S+(r,R,t+T)n„(r,R,t)]. 



70 



(6.24) 



From the random source method, we can evaluate ^m"(s, t+T)M"(s, t+T) from 



^«'^(s,t+TK(s,t+r)) = ^M-;(s,t+T;s,t+r) 

{a la 



(6.25) 



Since the configurations are gauge fixed to the temporal gauge before we estimate 
the correlator, the time-link link variables are unit matrices. Hence, we measure the 
product of the parallel transporters in Eq. (6.24) (i.e., Tr S^(i", t+T)Ilu[r, R, t) ) 
by two bent Wilson loops shown in Fig. 6.3. In terms of 



we have 



s{s,t+T) = ^M-;(s,t+r;s,t+r) 

a 

M(r,R,t,T) = 'Iv[E+(r,R,t+r)n„(r,R,i) 



Chab{T,R) = ^^^5(s,t+T)^/;(s-r)^(r,R,i,r). 



(6.26) 
(6.27) 

(6.28) 



We can express the wave function of the tttt state in momentum space, 

w{s, r) = w{s-r) = ^ e^P-^"-"V(p). (6.29) 
p 

Then we can rewrite Eq. (6.28) as 

Chab{T, R-) = ^ E [^(p) E t+T)e-'^-^ ^ ^(r, R, t, T)e'^-^^. (6.30) 
If we perform the following Fourier transforms. 



M(-p,R,t,T) = ^M{v,R,t,T)e'^-^ 



(6.31) 



Figure 6.3. An illustration of the calculation of JE. 
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S{p,t+T) = E'5(s,^+T)e-'P•^ (6.32) 

s 

we arrive at 

Chab{T,R) ^y^Y.T.^iP)S{p,t+T)M{-p,R,t,T). (6.33) 

Now we discuss the symmetry of the wave function w{s — r). It depends on 
the meson, which was created in the decay. For scalar meson production (with 
J^*^ — O"*""*"), the ^(s — r) must be in a n„ representation. This can be achieved by 
making w{s — r) odd in the x direction and even in the y and z directions, where x 
is the direction of the transverse extent of the Uu state. Another method to check 
this symmetry is to consider the parity, since Px, Py and CPz are conserved in the 
heavy decay in heavy quark hmit. Now we consider our n„ potential state, which 
is odd under P^ and even under Py and CPz- The operator is even under all 
three operations. Therefore, the symmetry requires that the wave function w be 
odd under the inversion P^ and even under Py and P^ (since C = +1 for the scalar 
meson). So, we choose the wave function as 

w{\i, r) = sin(A;a; • x) cos{ky ■ y) cos{kz ■ z), (6.34) 

where the k = {k^, ky, kz) is the chosen momentum. It have the desired symmetry. 
In terms of plane waves 

' r 

= ^|<^P,k+(^p,ki+(^p,k2+<^p,k3 - (^p -k-(^p ,-kl-(^p -k2-(^p -k3 |, (6.35) 

where ki = (A;^, ky, -k^), ka = {k^^, -ky, kz), and kg = {k^, -ky, -kz). 
We can then rewrite Eq. (6.33) as 

Chab{T, R)k = ^ E ^{siK t+T)Mi-k, R, t, T) - Si-k, t+T)M{k, R, t, T) 
+ 5(ki, t+T)M{-k^, R, t, T) - 5(-ki, t+T)iE(ki, R, t, T) 
+ 5(k2, i+T)iE(-k2, R, t, T) - S{-k2, t+T)M{k2, R, t, T) 
+ S{k^,t+T)M{-k^,Ii,t,T) - S{-k^,t+T)M{k^,Ii,t,T)^. (6.36) 

We use this formula to calculate the desired correlator Chab- 



CHAPTER 7 



HYBRID MESON DECAY FROM THE 



LATTICE 



For notational simplicity and without loss of generality, in this chapter we 
consider the generic strong decay 



where A and B are stable particles. 

In Chapter 5, we studied the hybrid meson decay channel. It is possible to 
employ some rather restricted conditions on a lattice to extract the transition matrix 
elements (x) [4]. In this chapter, we first use the time-dependent perturbation 
theory to study the Chab correlator. Then we introduce a method to extract x 
directly from lattice [4]. 

7.1 Perturbation expansion for Chab correlator 

In this section we use perturbation theory to study the correlator Chab in 
Sec. 6.3, that is to calculate the correlator from the 11^, state (i.e., H state) to 
two-body states (namely, AB state). To use this approach, we should discuss 
whether the perturbation theory here is justified or not. The condition we need is 
that r <^ A£', where F is the decay rate, and A£' = Eh — Ea- This is the narrow 
width approximation. In Chapter 8, we will discuss this question. 

For notational simplicity we refer to the lattice operator at the source H as 
(Please see details in Chapter 6), 



H^A + B, 



(7.1) 



OH{r,R,t), 



(7.2) 
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and we refer to the lattice operator at the sink AB as 

OABir,R,t+T) = Y,OA{r,R,t+T)wis-r)OB{s,t+T), (7.3) 

s 

where w{s—r) is the distribution wave function of the meson B. 
In this notation the time shce correlator is given by 

CHAB{t,R) = ^E(^AB(r,R,t+T)0^(r,R,t)), (7.4) 

where R indicates that the correlator CHAsit) is calculated at a fixed orientation 
of the QQ, and r is the midpoint of the QQ (see Fig. 6.2). The system can be 
described by a Hamiltonian H, that is, 

H^Ho + Vi, (7.5) 

where Hq is the unperturbed Hamiltonian, and Vj is the transition potential. In 
perturbation theory we assume that Vj is numerically small. Hence, 



^-Hot ^ J2 e-^<'(*-*iV,(t)e-^<'*i + • • • . (7.6) 

ti=0 



Therefore, we obtain 

Chab{T, R) = h E (o |e^(*+^)Oi5(r, R, 0)e-^(*+^) e-^'Onir, R, 0)e-^* 



^ r,t 



^ E E (O ^is(r, R, 0) e-^''^^-*^) Viih) e"^"*^ 0^(r, R, 0) 0) + 



T,t ti=0 



^ E (O R, 0) e-^»^ OH(r, R, 0)| o) (7.7) 

The second term in Eq. (7.7) is zero, because the Hq does not, by itself, cause the 
transition. Thus we can rewrite Eq. (7.7) as 



Chab{T, R) = ^ E E (o O^sir, R, 0) e'^^^^"*^) Vi{h) e"^"*^ 0^(r, R, 0) 



r,t ti=Q 

(7.8) 

In this dissertation, we model the transition potential as 

Vj{t') ^ E E E ^(R') (R'' r', ^') Cn„ (R', r', t') a^VV'(s', t') ^/(s' - r') , (7.9) 

R' r' s' 

where ,t(R') is the overall normalization factor that controls the transition ampli- 
tude, C!,+ is a creation operator for the E+ state, Cn„ is an annihilation operator 
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for n„ state, and wj{s'—r') is the distribution wave function of the ipi/j{s' ,t'), which 
is dictated by the decay. 



T-l 



CHABiT,R) = a;(R)-5:^i:i:i:(0Ol^(r,R,t')e-^°(^-*')Ct4R',r',t') 

T,t t'=0 R' r' s' ^ " 

Cn„(R', r', 0)V/iKs', f ) «'/(s' - r') e"^"*' e»^(r, R, 0) 



0). (7.10) 



If we plug into the operator Ol^g{r, R, 0), we obtain 



T-l 



Chab{T,R) = xiR)-J2Il E E(O^E+(r'R-'0)V^V'(S'^+^K(s-r)e 

" r,t t'=OR',r' s,s' ^ " 



-Ho(T-t') 



r', 0)CuM', r , 0)^^(s', t') wj{s'-r') e"^"* On„(r, R, 0) 



r,t t'=0 s,s' P,q 







R)^ E E E E ^« e" e-^n..',.c,.(s-s')^*(q)^^(q) 







* q t t'=0 s,s' 



V^^/;(s,t+T)V^V(s ,t +t) 0), 



(7.11) 



where is defined in Eq. (4.6), is defined in Eq. (4.12), and we express the 
wave function of the tttt state in momentum space. 



w(s-r) =Ee'P'^'"'^w(p) 
p 



(7.12) 



We can rewrite Eq. (7.11) as 

T— 1 

Chab{t,r) = x(R)y E ^*(q)^/(q) E ^ae"''^^ ^^"*'^Z,e-^""*' X 

t'=Q 

O) , (7.13) 



^ E E e^'*'^'"''^ (o |V'V'(s, i+r)V'V(s', + 1) 

t s,s' 

where V = N^NyN^. In Chapter 6, we calculate Cbb(T, k) through 

CBsiT, k) = ^ E E ^'''■^'"''^ (O t+T)V;V'(s', t) I O) . (7.14) 



t s,s' 



Therefore, we obtain 



T-l 



c^^s(r,R) = x(R)i^E^*(q)^Kq) E 



t'=0 



(7.15) 
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From Eq. (4.6) and Eq. (4.12), we have 

CHH{R,t') = e-^n-(«)^*'. (7.17) 

So we can rewrite this result in terms of all the correlators as 

CHAB{T,R)=J2w*{cDwi{q) J2 CAA{n,T-t')x{R)CHH{R,t')CBB{q,T-t'). 

(7.18) 

In Chapter 6, we discussed the symmetry of the wave function w, and we chose the 
wave function w as Eq. (6.34). Then we obtain 

Chab{T, R, k) = ^^^i^ X? C!aa(R, T-t')x{R)C„H(R, t')CBB{T-t', k), (7.19) 

where we consider the cubic symmetry (namely, Cbb(T— i', R, — k) = Cbb{T — 
i',R, k), etc). Here we absorb the i into ^/(k) (i.e., iwi{\<i) — >■ wi{\<i)). Note that 
the wave function wj must have the same symmetry as the wave function w. 

7.2 Decays from the lattice 

In practice it is possible to extract the transition amplitude x directly from the 
lattice [4]. Here we describe this approach in principle. Using a suitable lattice 
operator to create the hybrid state H &X, t — Q and annihilate a two-body state 
with relative momenta k and — k at time t. In Sec. 7.1, we used perturbation 
theory to obtain the correlator Cn-ABit) from an H state with mass Eh and a 
two-body state (^4, B) with energy Ea, Eb, that is, ^ 

ChabU, R, k) = ^^^^^ Caa(R, t-t')x{\i) C^^h(R, t') CBB{T-t', k), (7.20) 

where the summation is over the intermediate value t', and t' is an integer. By 
obtaining Zh from the correlator Cn^Hit) from H ^ H, and Za from the correlator 
CA-A{t) from A ^ A, also measuring the correlator CB-Bif) from B ^ B, and the 



^This equation is same as Eq. (7.19), just change T into t. 
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correlator Cn-ABit) from H AB, we can extract the transition amplitude a;(R) 
in principle, that is, 



^ ^ CaaCR: t-t') ChhCR, t') CBB{t-t\ k) 
t'=0 



If we define 

x{k, R) = x(R) wi{k) V, (7.22) 

we obtain 

x(k, R) = — CHAB{t,R,k)Z^Z, _ ^^23) 

Y,CAA(R,t-t') CHH(R,t') CBBit-t',k) 

t'=0 

Therefore, from Eq. (7.23), we can extract the transition amplitude a;(k, R) from 
the lattice. 

In this dissertation we focus on the decay channel to tttt branch, that is, 

^^Xfe + TT + TT, (7.24) 

where tttt is in scalar channel. In the lattice simulation, due to the taste symmetry 
breaking, there exist the splittings between pions (tt) of various tastes [60] [10] [63]. 
We classify the pions into five types, namely, tts, tt^, vry, tt/, and ttt which corre- 
spond the pions with taste pseudoscalar, axial, vector, singlet, and tensor respec- 
tively. In order to calculate the decay channel to tttt branch correctly, we should 
count these five channels correctly. In Ref. [64], the authors predict the weights 
for each taste of the tttt intermediate state in bubble contribution channel. We list 
them in Table 7.1. 

In our model, the transition amplitude of the Goldstone h^ttq can be calculated 
in lab frame (see details in Appendix D) by the formula 

lab 



xpp{p) = (7r5(pi)7r5(p2), A Vr H) , (7.25) 

where V}'^^ is the interaction potential in lab frame, which is denoted in Eq. (D.6), 
and the subscript PP specifies the Goldstone ttsTTs channel, that is, 

xpp{pi,P2) = (7r5(pi)7r5(p2) VV'(P) 0)x(p)5(^)(Pj^-Pa-p) 
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Table 7.1. The weights for different tastes of the tttt intermediate states. The first 
column is the taste, the second column is its weight. 



Taste 


Weight 


I (singlet) 


-3/4 


V (vector) 


4/4 


T (tensor) 


6/4 


A (axial) 


i/i 


P (pseudoscalar) 


1/4 



&pp(Pi,P2) x{p) 5^^\Ph - Pa - p), 



(7.26) 



where x{p) is calculated by Eq. (D.27) (see details in Appendix D), p = pi + p2, 
and bpp(pi,p2) = ^7r5(pi)7r5(p2) '0'0(p) 0^ is a constant, which does not depend 
on R, calculated to be [64] 

' ~ 1 



^'PP(P1,P2) = 



R2 
-^0 



1 



a' V 4L3 ^2^^^(pi) ^2£;,,(P2)^ 



(7.27) 



where Bq is the coupling constant, £'^p(p) = y^p^ + M^^ and Mu^ is the mass of 
Goldstone pion. Please see details in Ref. [64] or Appendix D. 
Similarly, another channels can be also measured. That is. 



xaa(pi, P2) = &aa(pi, P2) x{p) S^^\Ph - Pa - p) 

a^vv(Pi,P2) = byY{pi,p2) x{p)6^^\Ph -Pa- p) 

2;ii(Pi,P2) = bu{pi,p2)x{p)6^^\PH -Pa-p) 

a:^TT(Pl,P2) = bTT{Pl:P2)x{p)6^^\PH -Pa-P): 



(7.28) 
(7.29) 
(7.30) 
(7.31) 



where the subscripts AA, VV, II, and TT stand for the axial, vector, singlet, tensor 
channels respectively. 

1 



^aa(Pi,P2) = 
Vy(Pi,P2) = 

^'//(Pl,P2) = 



«'V4L3y^^;;;(^y^^;;;(^ 
1 1 



V 4L3^2K^(pi) ^2E,,{P2) 
a^\4L^ ^2E,M V2^..(P2) 



(7.32) 
(7.33) 
(7.34) 
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Ott(Pi,P2) = 



V 4L3 ^2K.(pi) V2^..(P2)' 



(7.35) 



where i?7rt,(p), -^ttaIp)? -^ttjIp), and -EtttIp) ai'e denoted like ^'^^(P). 

The transition amphtude of the tttt channel can be calculated by the formula 

(tttt, a jy,"^"! h) = x^4p)6^''\Ph -Pa-p) (7.36) 

where Vj^^ is denoted by Eq. (D.6), and the subscript tttt specifies the tttt channel. 
After some algebra, we obtain 



(7.37) 



where 67r7r(pi,P2) = (j^'^ip) ■0'0(p) 0/ is a constant, which does not depend on R. 
If we consider Table. 7.1, we obtain 

bl^iPi, P2) = |&pp(pi, P2)+&aa(Pi, P2)+&yy(Pl, P2) -|&/7(Pi, P2) + |&|r(Pl, P2)- 

(7.38) 



In continuum limit, we can estimate 



07r7r(Pl,P2) = ^1 



(7.39) 



where i?7r(p) is the energy of pion with momentum p. In practice, we can rewrite 
Eq. (7.37) as 

1 



pE,{p^)pE,{p^y 



where 



a3 V 4L3 



(7.40) 



(7.41) 



We can also define a constant, that is. 



a3V4L3- 

In Chapter 8, we show how to calculate the decay rate from X7r7r(p). 



(7.42) 



CHAPTER 8 



OUR NUMERICAL RESULTS 



In this chapter we give our lattice numerical simulation results. First we present 
the calculation of the mass of the /o meson. Then we show the detailed procedure 
to obtain the lattice transition matrix element and calculate the decay rate. 

In this dissertation, we use 520 dynamical MILC fermion configurations with 
Nf — 2 + 1 flavors of sea quarks (namely, amu,d — 0.005, and am^ = 0.05). The 
coupling constant is ^ = 6.76. The lattice dimension is 24^64, with the lattice 
spacing a of around 0.12 fm. All the results we show in this chapter come from this 
ensemble of gauge configurations. 



Recently there exists a growing interest in the /o meson. The analyzes of the 
TTTT scattering phase shift suggested a brand resonance of the /o meson with 7 = 
and J^^ — 0++ [65]. In this section we first discuss an interpolating operator for 
the /o meson. Then we give the detailed procedures for determining its mass. 

8.1.1 Operator for /o meson 

We choose an operator with isospin 7 = and J^*^ = O"*""*" at the sink and source, 

- , . _ ^ Ur{x)Ur{x) + dr{x)dr{x) . . 



where r is the index of the taste replica, rir is the number of the taste replica, and 
for notational simplicity we omit the color index and spin index. 
The timeslice correlator for the /o meson is given by. 



8.1 /o meson 




(8.1) 




(8.2) 
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where the x — (x, t) is lattice position, is zero vector (namely, = (0, 0, 0)), (/o) 
is the vacuum expectation value, and k is the chosen momentum. 

In Chapter 6, we discussed this correlator in detail. Prom Ref. [60], we know 
that rir = 1/4 for the 1 + 1 + 1 theory. Thus, we obtained 

C(t,k) = C^^(t,k)- V^(t,k), (8.3) 

where C{t,k) is described by Eq. (6.6), C"^'^(i,k) is described by Eq. (6.8), and 
C^'^(t, k) is described by Eq. (6.7). The superscript CC stands for the connected 
contribution. The superscript DC stands for the disconnected contribution. 

8.1.2 Procedures for determining the /o mass 
As we discussed in Chapter 6 we used the point sources and the point sink 
operators for the /o propagator. For the point source, we just set a 1 on only 
one point on a chosen time slice, and zero everywhere else. In order to improve 
the statistics, when we calculate the connected part, we computed correlators from 
eight source time slices evenly spread through the lattice (i.e., only one source time 
shce was chosen at a time), and averaged the correlators. Moreover, for the nonzero 
momentum mesons we used a quark source with 1 and an antiquark source with e*'''^ 
on only one point on the chosen time slices, where the k is our chosen momentum, 
and the s is the position of the quark source (see Fig. 6.1). 

From Ref. [64] , we know, for staggered quarks, that the meson correlators have 
the generic single-particle contribution [38] 

C{t) = ^ A (e-'"'* + e-'"'^^*-*)) + ^^;(-l)* (e"'"'* + e-™*^^*-*)) , (8.4) 

i i 

where the oscillating terms correspond to a particle with opposite parity. In our 

case for the /q meson correlator, we assume only one mass with each parity in the 
fits of Eq. (8.4). ^ From the discussion in Appendix B, we must consider the two 
meson contributions. Therefore, we have 

C(t) = Ai (e-"^i* + e-"^^(^*-*)) + A2{-iy (e""^^* + g-^^^^^*"*)) + B^ff'^{t), (8.5) 

^In Ref. [64] we show that, in our concrete calculation, our operator is the state with 7(g) 7 and 
state of 7o75 (S) 7075 (its oscillating term). 
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where Bf^^^(t) is denoted in Eq. (D.14). If we rearrange the terms in Eq. (D.14), 



we obtain 



B 



SxPT 



it) 



fBubble(^) + 5vfvhairpin(^) + ^^fAhairpin(i) ) ■, 



(8.6) 



here Bq is the couphng constant, 5v — cl'^^'v the hairpin couphng of a pair of 
taste-vector mesons, 5a — o?5'j^ is the hairpin couphng of a pair of taste-axial 
mesons, and 



^Vhairpin 



k 



-4- 



(p+k)^+M2^ + ^k2+M^^)t 



-(^(p+k)2- 



Tj'V 



-(^(p+k) 



-Cy Cv, 



(p+k)2+M2 +^i?TM^ 



+ 



(8.7) 



iBubble 



(t) 



^ ^-(^^(p+k)2+M^2^+^/k5TiW|;)t ^-(y/(p+k)2+M2^+^k2+M2jt 



E 



16 



6=1 



9 y'(p+k)VM2,^kHM2, ^(p+k)^+M2^^k2+M; 

(p+k)2+M2^ + ^k2+M2jt 



2 



(p+k)2+M^,+Jk2+Afi )t 



+ 



^(p+k)VM2^^k2+M2^ ^(p+k)^+Mi^^k2+M| 



, (8.8) 



where Sy, Sa, Bq, Cy^, Cy^,, Ca^,, Ca^, are constants, which are given in Appendix 
B, M^y, Mr^Ai Myf'v, Mjj/A are given in Appendix B, and for fAhairpin(i), we just 
require V ^ Am fvhairpm(i)- 

The tree-level masses of the mesons are [60] 

M]f,^ = Boinif + ruf) + a'^Af,, (8.9) 

where //' arc two flavors which make up, b = 1, ... 16 arc the taste, and the term 
of a^Aft comes from taste symmetry breaking. Here for brevity, we denote 

(8.10) 



By convention the Mup is the mass of the Goldstone pion, and the Mkp is the mass 
of the Goldstone kaon. 
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All the masses needed for Eq. (8.7) are well known. From the mass-squared 
splittings in Table III of Ref. [66] also the aMup = 0.1594, aMxp = 0.36523 in 
Ref . [10] , we obtain 

aMu^ = 0.2342, aMu^ = 0.2694, aMu^ = 0.2966, aMu^ = 0.3205. 
uMk^ = 0.4036, aMKr = 0.4250, aM^^ = 0.4428, uMkj = 0.4591. 

The is estimated from the Gell-Mann-Okubo formula 

ml = (8-11) 

that is, 

aM^j = 0.4958. (8.12) 

The first term [J^'-^ — 1++) in Eq. (8.5) gives the mass of the /o meson, and 
the oscillating term [J^'" = 1"' ) in Eq. (8.5) gives the mass of the rjA meson. In 
Ref. [64], we do a fit to all of the correlator data at once. In this grand fit, the 
masses of the /o meson, oq meson, 7705 meson, and ttos meson, the couplings of Sy, 
6a, and Bo were determined. We choose all the data within the fitting range of 
t from 5 to 18 (namely, the minimum distances tmin — 5). Our fitting gives an 
acceptable result with a of 127 with 114 degrees of freedom, with 

am/o = 0.47(7) (8.13) 

Prom Ref. [10] , for our set of configurations, we know 

a"^ = 1.634(45) GeV, (8.14) 

then gets 

m/o = (768 ± 136) MeV. (8.15) 

Also we get 

a^o = 3.12 ± 0.42. (8.16) 
In next section we use this number to evaluate the decay rate. 
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8.2 Lattice transition matrix element 



Prom the discussion of the symmetry of the wave function in Sec. 6.3, we 

know the possible minimum momentum is k = (1,0,0)— for the scalar meson 

Ij 

creation [4]. In this section, we will show the detailed procedures to obtain the 
overall decay rate at this chosen momentum. Then we will consider the momentum 

of k = (1, 1, 0) — , k = (1, 1, 1)—, and k = (2, 0, 0)—. From the data of these 
four chosen momenta, we can obtain a decay rate that has physical meaning. 



In Chapter 4, we used the LBO approximation to solve for the gluonic field 
around a static quark- antiquark QQ at separation R to determine first the static 
potential. Then we solved the Schrodinger equation with that potential to deter- 
mine its spectrum and wave function. Prom Pig. 4.11, we note that the potential 
(n^,) that binds the hybrid meson is comparatively flat. Hence, the hybrid meson 
has a broad radial wave function, which we denote hj Uu va. Chapter 4. This 
observation allows a quick estimate for the decay rate. For example, in Chapter 5, 
we justified the decay from the hybrid meson state to the P-wave Xh state (i.e., we 
refer to the wave function of this state as u^^ in Chapter 4). In order to integrate 
over R easily, we rescale Pig. 4.11 in terms of lattice spacing a, and show it in 
Pig. 8.1. We see that the lack of the nodes in the relevant wave functions results 
in a quite large spatial wave function overlap factor. Moreover if we assume the 
transition rate is independent of i?, we obtain 



where uu^ u^^, are normalized. Prom the other side, this big overlap factor justifies 
the decay from the hybrid meson state to the P-wave Xh state. Later in this chapter, 
we calculate the decay rate of this channel. 

8.2.2 Simulation procedure 

In Chapter 6 we have already discussed how to create operators for the n„, 
-I- 5'(0"'""'"), and we use Eq. (6.33) to calculate the correlator Chab- We also 



8.2.1 The wave function of the excited gluonic states 




(8.17) 
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0.0 







2 



4 



R/a 



Figure 8.1. Reduced radial wave function versus R in units of a. The solid curve, 
dotted curve, and short dash curve specify the reduced radial wave functions for 
the 15, IP states of S+, and n„ hybrid state, respectively. 

discussed the symmetry of the wave function w{y). 

In practice we evaluate the difference of two Wilson loops. As we discussed in 
Chapter 6, we call this observable iE(r) and its spatial Fourier transform iE(k). 
The fermion loop of light quark meson (/o meson) is evaluated at each spatial point 
s at time t + T and its Fourier transform is 5'(k). We actually used 20 iterations of 
APE smearing for both spatial ends of M. The transverse extent of the n„ end is 
two lattice spacings. We summed over all possible orientation to improve statistics. 

8.2.3 Our numerical results for decay rate 

For notational simplicity, we focus on the chosen momentum k = (1,0,0)—. 
For scalar decays, we choose the wave function as in Eq. (6.34). According to 
Eq. (6.36), in terms of these components of the momentum, we estimate the 
required correlator CnABit, R, k). From Eq. (7.23), we evaluate the bare transition 
amplitude x(k, R, t) at each t and R for fixed momentum k, that is. 



x(k,R, t) 



CHAsit, R, k) Za Zh 



(8.18) 



Y,CAA{R,t-t') CHHiR,t') CBB{t-t',k) 



t'=0 
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We measure from the correlator Cn-nit) from H ^ H a.s explained in 
Sec. 4.4.2, and Za from the correlator CA-A{t) from A ^ A a.s explained in 
Sec. 4.4.1. The Za and Z^ values for different R are shown in Fig. 8.2, where 
the burst one is for Z^, and the octagon one is for Za- 

When we evaluate the bare transition amplitude a;(k, R, t), we use jackknife 

method. Here we should point out when we construct the jackknife sample, we just 

throw out one measurement. 

For illustration, the values of a;(k, R, t) for — = 1 and momentum k = (1, 0, 0) — 

a L 

at different t are shown in Fig. 8.3. For the range from t = 2 to t = 5, there exists 
a plateau for x(k, R, t). We fit this plateau to get x(k, R) = 0.3253 ± 0.0093 with 
a of 0.03 for 3 degrees of freedom. 

In our lattice simulation we obtain the x(k, R, t) values shown in Fig. 8.4 and 

Fig. 8.5. From these data we can extract the values of x(k, R) for momentum 

271 2tt 
k = (1,0,0)— at all separations R. In Table 8.1, for k = (1,0,0)— we list all 
Ij ij 

the values of x(k, R) for different i?, fitting range, etc. Also we plot the values 
of x(k, R) in Fig. 8.6. From Fig. 8.6, we note an increase of x(k, R) with 



1.5 



1,0 



ISI 



0.5 




0.0 I ' ' ' ' ^ ' ' ' ' ^ ' ' 1 

5 10 

R/a 

Figure 8.2. The overlap amplitudes of Za and Zh versus i?, the burst one is for 
Zh, and the octagon one is for Za- 
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Figure 8.3. The x(k, R, t) values for — = 1 versus t with k = (1, 0, 0)^. 
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Figure 8.4. The x(k, R, values with momentum k = (1,0,0)— versus t. Here 

Ij 

R 

values for — = 1, 2, ■ ■ ■ , 6 are represented by the plus sign, cross, diamond, square, 
a 

octagon, and fancy diamond, respectively. 
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27r 

Figure 8.5. The x(k, R, values with momentum k = (1,0,0)— versus t. Here 

Ij 

R 

values for — = 7, 8, ■ ■ ■ , 12 are represented by fancy square, fancy cross, fancy plus, 
a 

burst, cross, and dollar sign, respectively. 



Table 8.1. The x(k, R) values with momentum k = (1, 0, 0)—. The first column is 

the separate distance, and Column two is the transition amplitude. The remaining 
columns are the time range for the chosen fit, and number of degrees of freedom 
for the fit. 



R/a 


ax 


Fitting 


■ range 


XVD 


1 


0.3215(105) 


2- 


-5 


0.0266/3 


2 


0.5430(196) 


2- 


-5 


1.81/3 


3 


0.7170(272) 


2- 


-5 


2.51/3 


4 


0.8636(303) 


2- 


-3 


3.71/1 


5 


0.9520(382) 


2- 


-4 


2.02/2 


6 


1.1328(826) 


3 


-5 


3.01/2 


7 


1.270(119) 


3- 


-5 


0.987/2 


8 


1.511(163) 


3- 


-5 


1.07/2 


9 


1.615(189)) 


3- 


-5 


1.36/2 


10 


1.567(210) 


3 


5 


3.51/2 


11 


1.672(297) 


3 


-5 


0.893/2 


12 


1.465(395) 


3- 


-5 


0.16/2 
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Figure 8.6. The x(k, R) values with momentum k = (1, 0, 0)— versus R. 

Ij 

respect to R. and a hint of saturation at large R. One approach to interpret the 
i?-dependence of a;(k, R) is by noting that the scalar meson wave function w has a 
node at the center of the 11^ state in the transverse direction, since the w is odd 
in this direction, and therefore it is sensitive to the transverse width of the excited 
gluonic flux [67] in the n„ state. Just as we find, this can explain the increase [67] 
with longitudinal extent R leading to a plateau. 

Now we have the bare transition amplitude values a;(k, R). These values depend 
on R. In Appendix D, we develop a formula to integrate out R. From Eq. (D.35), 
we obtain the overall transition amplitude for fixed momentum k (i.e., 5(k)), 



We obtained the values of uh{R) and Uy^^^{R) from the LBO treatment in Chapter 
4, and we also plot them in Fig. 8.1. From Fig. 8.1, we note that when R > 3a, 
the value of u-^^{R) is approximate close to zero. Hence, only values of x(it!, k) for 
R — la, 2a, and 3a mainly contribute to the x(k) value. 




(8.19) 
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27r 

Table 8.2. The x(k, R) values with momentum k = (2, 0, 0)—. The first column 

Ij 

is the separate distance, and Column two is transition amplitude. The remaining 
columns are the fitting range, and number of degrees of freedom for the fit. 



R/a 


ax 


Fitting range 




1 


0.4007(114) 


2-5 


0.952/3 


2 


0.6433(181) 


2-5 


1.8/3 


3 


0.8515(221) 


2-5 


0.658/3 


4 


1.0001(242) 


2-5 


3.85/3 


5 


1.0006(276) 


2-4 


1.51/2 


6 


1.2160(310) 


2-4 


2.52/2 



Prom the values of x(i?, k), UH{R),u^i,iR), we can get an overall x value for 
27r 

momentum k = (1, 0, 0)— by Eq. (8.19), that is. 



ax(lOO) = 0.1693(60), 
where the 100 specifies the chosen momentum. 



(8.20) 



From the same procedure, we can get overall x value for momentum k = 

271 

(2,0, 0)— , that is, 

Ij 



ax{200) = 0.2007(50), 

2% 

the overall x value for momentum k = (1, 1, 0) — , that is, 

Ij 

ax(llO) = 0.1482(41), 
27r 

and the overall x value for momentum k = (1, 1, 1) — , that is, 

Ij 

ax{ni) =0.1115(28). 



(8.21) 



.22) 



(8.23) 



For concreteness, in Tables 8.2, 8.3, and 8.4, we list the values of x(k, R) for 

2% 27r 27r 

different it! with momentum k = (2, 0, 0) — , k = (1, 1, 0) — , and k = (1, 1, 1) — 

Ij Ij Ij 

respectively. Since when R > 3a, the values of u^^{R) are negligible, and only 
x{R, k) of -R = la, 2a, 3a mainly contribute to the a;(k) value. We omitted the 
results for > 6a in the table. 
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Table 8.3. The a;(k, R) values with momentum k = (1, 1,0)—. The first column 

is the separate distance, and Column two is transition amplitude. The remaining 
columns are the fitting range, and number of degrees of freedom for the fit. 



R/a 


ax 


Fitting range 




1 


0.2909(85) 


2-5 


0.344/3 


2 


0.4695(150) 


2-5 


0.841/3 


3 


0.6352(174) 


2-5 


3.1 /3 


4 


0.7512(217) 


2-4 


3.35/2 


5 


0.9205(483) 


3-5 


0.0467/2 


6 


1.0056(620) 


3-5 


0.0896/2 



Prom our grand fit in Ref. [64], we found 

aSo = 3.12 ±0.42. (8.24) 
Hence, from Eq. (7.38), we have 

= 0.0230(31). (8.25) 

Therefore, from Eq. (D.37), we obtain the transition amplitute of the tttt channel 

2n 

for momentum k = (1, 0, 0)—, that is, 

ij 

a^x{im) = 0.00389(66), (8.26) 

where the 100 specifies the chosen momentum n = (1,0,0). Prom the same 
procedure, we can get overall transition amplitude x for other three momentum. 
We list all the results of transition amplitude in Table 8.5, where the first column 
is the chosen momentum, and the second column is overall transition amplitude. 

Prom these x-values, in Appendix D we relate the overall LBO decay amplitude 
x{p) to the lattice transition amplitudes x.^T^{\i) as Eq. (D.38). That is, we can use 
the follow formula 

^7r7r(p) = -4 H H o?Ns3i{pr)Yii{e, 0) x(k) sin(/c^x) cos{kyy) cos{k,z). 

(8.27) 



91 



Table 8.4. The x(k, R) values with momentum k = (1, 1, The first column is 

the separate distance, and the second column transition amplitude. The remaining 
columns are the fitting range, and number of degrees of freedom for the fit. 



R/a 


ax 


Fitting range 




1 


0.2742(65) 


2-5 


0.0458/3 


2 


0.3843(105) 


2-5 


0.672/3 


3 


0.4493(117) 


2-5 


1.85/3 


4 


0.3425(132) 


2-4 


5.11/2 


5 


0.3540(308) 


3-5 


2.41/2 


6 


-0.00207(436) 


3-5 


2.79/2 



Table 8.5. The overall transition amplitude a;. The first column is the chosen 
momentum, and Column two is overall transition amplitude. 



k 


a^x 


(l,0,0)27r/L 


0.00389(66) 


(2,0, 0)27r/L 


0.00461(74) 


(l,l,0)27r/L 


0.00341(55) 


(l,l,l)27r/L 


0.00256(41) 



This result gives the transition amplitude x at any given momentum p. Here Ng 
is a normalization factor, ji{pr) is the spherical Bessel function, and Yi ms(^,0) is 
the spherical harmonic function. In terms of our measured values, we have 

XnAp) = a^Nsji{pr)Yu{9, 0) x |£(100) sin(a;) + f (200) sin(2a;) 

-|- £(110) sin(a;) cos(|/) + £(111) sin(a;) cos(|/) cos(2;)|, (8.28) 

where we sum over just a few small values of k. With these values in Table 8.5, 
we evaluate Xj^^ip) at any given momentum p. We show the transition amplitude 
Xj^T^{p) versus momentum p in Fig. 8.7. Since we just measure a few x's with small 
momentum. As discussed in Appendix D, the values shown in Fig. 8.7 for high 
momentum are not reliable. We test that it is reliable for the small momentum p. 
For the decay 

^^Xb + TT + TT, (8.29) 
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0.010 



0.008 



0.006 



XI 0.004 




Figure 8.7. The transition amplitude x{p) versus momentum p. 

Prom Fermi's Golden Rule for the transition rate F, we need to calculate the phase 
factor. Because the calculation of phase factor is tedious. We show the detailed 
procedure in Appendix E to obtain a formula to cvahiatc the transition rate F, 



r = Stt^ (^^) J^^^^ dp p" x^^ipf 



^.30) 



^4 (AM)2-p2' 

where Ma is the mass of the Xb state, Mh is the mass of the H state, AM = 
Mh- Ma, and 

p^ax = ^J{AMy - Ami. (8-31) 

Here 171^^ is the mass of pion meson. 

For our simulation, we have Mh = 10.828 GeV and Ma = 9.877 GeV, hence 
AM = 0.951 GeV, thus p^ax = 0.83 GeV. We can estimate x^^{p) by Eq. (8.28). 
Hence, by integrating Eq. (8.30), we can obtain overall transition rate 



F = 3.62(98) MeV, 



(8.32) 



for the decay — > + tt + tt. In Chapter 7, we use the perturbation expansion for 
Chab correlator to find a way to extract the transition amplitude x. To use this 
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method, wc should discuss whether it is justified or not. The condition we need is 
r <^ A£^. Since F <^ A£^, our lattice method for decay rate is justified. 



CHAPTER 9 



CONCLUSION 



In this chapter, we first discuss the evidence for the /o meson. Then we discuss 
the evidence for Ught quark hybrid exotic mesons. Next we discuss hybrid exotic 
quarkonium decay. Finally, we summarize our results. 

9.1 Evidence for /o meson 

For the past 45 years, the question of the existence and nature of the /o-meson 
has attracted the attention of many authors [68] . The subjects of investigation are 
the internal properties of the /o-meson and its role as an intermediate particle in 
various processes, both in vacuum and in hot and dense matter. It appears as a 
broad enhancement in the isospin / = 0, J = tttt scattering amplitude. The 
experimental value of the /o-meson mass is not accurately determined and lies in a 
wide interval [12] 

M/, = 400 - 1200 MeV, (9.1) 

This uncertainty is explained by large values of the decay width of this meson into 
two pions [12] 

T/^ = 600- 1000 MeV. (9.2) 

Recently, the Fermilab experiment E791 [69] showed strong evidence for the 
existence of an / = 0, J = 7r"'"7r~ resonance in charmed meson decay. 
The 7r"'"7r~ resonance in E791 is compatible with the scalar meson /o(500). It was 
also observed in the Cabbibo-suppressed decay — > 7r~7r+7r+ [69]. The D'^ — > 
/o (500)71+ decay contribution is dominant, being responsible for approximately half 
of the decays thought the particular resonant sequence: — >■ tt -|- /o — )■ 7r~7r+7r+ 
decay [69]. The measured mass and width of the /o meson are 478 ± 31MeV and 
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338 ± 48 MeV, respectively. Several other experimental results can be explained by 
the existence of the /o meson resonance [70] [71] [72] [73] [74]. 

9.2 Evidence for light-quark hybrid exotic mesons 

There are two ways to study hybrid exotic mesons. One approach is to use lattice 
QCD calculation. A summary of current LQCD calculations for the lowest lying 
hght quark {u or d) J^*^ = 1~+ hybrid exotic meson is hsted in Table 9.1. These 
lattice QCD calculations are quenched or partially quenched. All of this numbers 
indicate that the rough mass of hybrid exotic meson is about 1.8-2.1 GeV. 

Another approach is to directly perform experiment [79]. Several high statistics 
experimental searches for hybrid exotic mesons have been carried out. These 
experiments show credible evidence for the existence of J^'-^ = 1~+ hybrid exotics 
7ri(1400), 7ri(1600), 7ri(2000). The 7ri(1400) was first reported by the VES Collab- 
oration at IHEP [80] and was confirmed by the E852 Collaboration at BNL in the 
7r~p — )■ rjTi'p reaction [81]. This was followed with additional confirmation by the 
Crystal Barrel Collaboration at the CERN LEAR facility in antiproton-neutron 
annihilation [82]. In the Particle Data Group (PDG) hstings [83], the mass of 
the 7ri(1400) is M = 1376 ± 17 MeV and its width F = 300 ± 40 MeV with 
observed decays to nrj. A second J^*^ — 1~+ exotic hybrid candidate, the 7ri(1600), 
was first observed at BNL E852 in 7r~p pp ^ 7i~^7i~7i~p [84]. Since that time 
further evidence for the 7ri(1600) has been provided through its decays into rj'ii [85], 
/iTT [86], and bin [87]. The PDG [83] lists the mass and width of this state as 
M = 1596i?t MeV and P = 312^^^ MeV. 

Table 9.1. Quenched lattice QCD calculations for the mass of uug hybrid exotic 
meson {J^^ — 1'"^). 



Collaboration (Year) 


Computed Mass (GeV) 


Reference 


UKQCD (1997) 


1.87 ±0.2 


[75] 


SESAM (1998) 


1.90 ±0.2 


[7(3] 


Mei and Luo(2003) 


2.013 ±0.026 


[77] 


MILC(2004) 


1.792 ±0.139 


[78] 
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Another candidate J^*^ = 1 + meson is 7ri(2000), which has been seen through 
its decays into bin [87] and /itt [86] at BNL E852. The J^^ for bi is 1+", and the 
J^^ for /i is 1"^+. This state is the least controversial of the announced exotics, 
because its mass is consistent with theoretical expectations and its decay modes are 
consistent with expectations for hybrid mesons. However the statistical accuracy 
of the existing data still needs improvement. 

Because the lattice simulation of the heavy quark hybrid exotics (namely, bbg) is 
difficult, and searching them from experiment need pretty high energy, hence there 
are no papers about this kind of heavy quark hybrid exotics. The result of this 
dissertation will bring a little contribution in this field. 

9.3 Hybrid exotic quarkonium decay 

Recent observations of charmonium states in S-meson decays suggest that char- 
monium hybrid mesons (he) ^ can be produced in B meson decay [89, 90, 88]. Some 
of these states are likely to be narrow with subsequent decays to J/ipn^n~ , etc. 

Recent developments in both theory and experiment lead us to expect that 
charmonium hybrids will be produced in B decays. The partial widths for B — > 
cc + X, with cc representing specific final states such as J/V', Xco, Xci, Xc2, etc. 
Possible decay schemes are [88]: /ic(0"'"~, 2+~) — > J/i}+ {7r~^ ,r},ri'), hc{l~~^) — >■ 
Vc + (7r'''7r~, 77, 77'). These decay modes can give distinctive signals. Both 0"*"" and 
2"*"" should decay via the J /ipTiTi cascade [88]. The hcil'^) state is expected to be 
the lightest exotic cc hybrid, and therefore in this case the cascade goes to jycTTTr. 

The X(3872) is a narrow state decaying into n'^n'J/ip, with a mass Mx ~ 
3872 MeV. Given the observed mass, in the charmonium region, it is natural 
to assert that the X(3872) is itself a charmonium state. However it is difficult 
to identify the X(3872) with any of the expected narrow cc mesons, leading to 
suggestions that it is a more exotic particle [91] [92] [93]. 



The notation he is introduced in Ref. [88]. I use his notation in the discussion of this section. 
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9.4 Summary 

We reviewed arguments that, in the heavy quark Umit, the decay channel of a 
J^^ — spin-exotic hybrid meson will be dominated by the creation of a flavor 
singlet light quark-antiquark. In this dissertation, we explored this decay in lattice 
QCD with Nf — 3 flavors of dynamical sea quark with masses close to the physical 
values. We studied the transition between an excited hybrid gluonic state (i.e., 
IIu state) and a ground state gluonic state (i.e., Xb state) with the production of 
a flavor-singlet light quark-antiquark pair. The dominant decay channel should 
be the tttt channel. We introduced a phenomenological model for the decay and 
directly studied the transition to two pions on the lattice. We obtained a little bit 
smaller transition amplitude for this decay channel. 

Our numerical result for the scalar /q meson is am = 0.47(7) (or mj,) = (768 ± 
136) MeV for this lattice spacing a = 0.12 fm, the mass of up quark {arriu = 0.005), 
and the mass of strange quark (am^ — 0.05). The energy release from the hybrid 
meson at 10.833 GeV to the Xb state at 9.893 GeV is 0.940 GeV. Therefore, it 
is impossible to decay to a scalar meson heavier than 940 MeV. For b quarks the 
relevant transitions is if —>x;, tttt. Our raw lattice result is an overall decay rate of 
3.62(98) MeV for the scalar meson. 

The evaluation of this transitions amplitude in this dissertation is a fascinating 
step for the phenomenological study of hybrid decays. Our results are consistent 
with our expectation that the spin-exotic hybrid meson states are relatively narrow. 
We hope our results will be able to guide the experimental search for a bbg and /o 
meson, etc. 

In this dissertation, wc performed our research on one ensemble of lattice gauge 
configuration. That means our results arc at a single lattice spacing and light quark 
mass. As we discussed in Chapter 3. The physical one should be in continuum limit 
(lattice spacing a — > 0). Hence, we still work on other different configurations. In 
the continuum limit, we get the decay rates that have physical means. 



APPENDIX A 



THE PC OF THE HYBRID MESON 

In this appendix, we want to prove Eq. (4.34), and Eq. (4.35), which we use in 
Chapter 4 to evaluate the PC of the hybrid meson. 

Here we first introduce the azimuthal angle 9. the zenith angle ijj of the axis 
of the QQ and the distance (radius) between Q and Q with respect to a fixed 
system of coordinates ^,r]X- Besides the fixed system of the coordinates ^, 77, C) 
we also introduce a moving system of the coordinates x, y, z. And the two system 
have the same origin. The z-axis is in the direction of the axis of QQ, and the 
X-axis is lying in the ^ — 1] plane. In bpdy-fixed system of the coordinates x, y, z, 
we also introduce the cylinder coordinates: radial p, angle a, and z. The angle a 
is measured in xy-plane from the x-axis with < a < 27r. 

The coordinates Cg^Vg^ Cg the gluon and the coordinates Xg, yg, Zg of the gluon 
have following relationship 

Xg — —^g sin ip + r]g cos ip (A.l) 
yg — —^g cos 9 cos t/^ + rjg cos ^ sin -0 + Q sin 9 (A. 2) 

Zg — —^g sin 9 COS i/j + % sin ^ sin + Cg cos 9. (A. 3) 

We know, under the operation of parity, that the 9, ip change according to 

9 ^ 7r-9 (A.4) 

tl) TT + t/;. (A. 5) 

Also the ^,rjX change sign. Prom Eqs. (B.1,B.2,B.3), we obtain that, under the 
operation of parity, the x and z do not change sign, and y change sign. This means 
that, under the operation of parity, 

a — >■ —a. (A. 6) 
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The wave function of the hybrid meson is, then, 



Xn,A,L,ML = 0n,A,i(P'«'^)^n,AAMi(r)6L,Mi,A(6')$Mi(^), 



(A.7) 



where L is the total angular momentum of the hybrid meson, M^, is the ^-component 
of L, A is the z-component of the angular momentum of the gluon, and n stands 
for the others quantum numbers, which specify the state of the hybrid meson, 
(f>n,A,L{p,C(, z) is the gluon wave function, and S^ a.l.MlI^") is the "radial part" of 
QQ wave function. The 9 dependence is given by OlMlaI^)- Finally, ^M^iip) gives 
the 'V dependence. 

According to the solution in Page 299 of Ref. [53], we obtain 



e 



(0) = i-i) 



(2L+1)!(L+Mi) 



cos 



X 



\ (L-A)!(L+A)!(L-M^)! (i+^^s 9) 

r. \ L-Ml 

^ (1-cos ^)^-^(l+cos ^)^+^ 



2 



X 



[ \dcos 9 J 

which we are called the left-handed gluon state, and 



(A.8) 



{ ( d 



(2L+1)!(L+Ml)! 



'1— cos 



\ {L-K)\{L+K)\{L-Ml)\ (1+cos 9)- 

\ L-Ml >> 



X 



:i-cos 



cos 



vL-A 



(A.9) 



[ \dcos 9 J 

which we are called the right-handed gluon state. It is obvious that the left-handed 
gluon state and the right-handed state are not eigenstates of parity. Now if we 
combine them, we can construct eigenstates of parity. 



{2L + 1)\{L+Ml)\ 
\ {L-A)\{L+A)\{L - Ml)\ 



X 



D 



L-Ml 



(1+cos 9)^ V^cos 9^ 
(1— COS 9) 2 id] 



1-cos ^)^(l+cos + 



(1-COS e)^-^( 1+COS 9) 



(1+COS 

where e = ±1. 



1^ lacos e 



L+A 



(A.IO) 
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Hence, we obtain 

P {©l,m^,a(^) + eeL,ML,-Am = ^L,M„A{n -e) + sOlmlM^ - (A.ll) 
Now if we consider the identity cos(7r — 9) = — cos 9, we obtain 

p{eL,M,AO) + ee^,M.MO)} = (-0^ 



{2L+1)\{L+M,)\ ._^^L-M, 

^ (L-A)!(L+A)!(L-Mi)! ^ ' 



1 cos^) ^ '''(l+cose)^-^(l-cose)^+^ ^. (A.12) 



(1+cos 9)^^ \9cos 9 J 
Hence, we obtain 

(A.13) 

The part of the gluon wave function (namely, $(^)) is 

^M.W^-^e^'"-^. (A.14) 



It is obvious that 

P^MAi^) = {-1)'''^mA^)- (A.15) 
And the "a" part of the gluon wave function in 0n,A,L(p, ck, z) is 

0A(a) = 4=e^^^ (A.16) 



It is obvious that 

P(/)A(Q;) = 0_A(a). (A.17) 

If we consider the intrinsic parity of the meson is -1 [14], the total parity of the 
hybrid meson is given by P = e{—l)^~^^ x (—1)^^ x (— 1)S that is, 

p = e(-l)^+^ (A.18) 

Let us define the PC quantum number of the gluon to be r], if we also con- 
sider that the symmetry of the spin wave function under charge-conjugation is 
(—1)'^"'"^ [14], we obtain the charge-conjugation of the hybrid meson, namely 

C = er]{-lf+^. (A.19) 



APPENDIX B 

BUBBLE CONTRIBUTION FOR ao 
CORRELATOR 

The mass of the oq meson can be rehably determined on the lattice. In order to 
determine the mass of the ao meson, we evaluate the ao correlator 

C{t) = 5^((i(x,t)M(x,t)M(0,0)rf(0,0)) . (B.l) 

X 

The extraction of the mass of the ao meson ( = 0+ and / = 1) is straightforward. 
However there exist many multihadron states with = 0+ and 1 = 1 which can 
propagate between the source and the sink. Of special interest in multihadron 
states is the intermediate state with two pseudoscalars P1P2 which we refer to as 
the bubble contribution (B) [63]. If the masses of Pi and P2 are small, the bubble 
contribution B gives a considerable contribution to the ao correlator, and it should 
be included in the fit of the lattice correlator in Eq. (B.l), that is, 

C{t) = Ae-'^''* + B{t) , (B.2) 

where we omit the unimportant contributions from the excited ao meson and other 
high order multihadron intermediate states. 

B.l Coupling of a scalar current to pseudoscalar 

Before we embark on the calculation of the bubble contribution, in this section 
we first derive the coupling of a point scalar current dr{x)ur{x) to a pair of the 
pseudoscalar fields at the lowest energy order of the staggered chiral perturbation 
theory (S^PT), where the subscript r in the expression Ur{x) is the index of 
the taste replica for a given quark flavor u. The effective scalar current can be 
determined from the dependence of the lattice QCD Lagrangian and the staggered 
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chiral Lagrangian on the spurion field Ai, where Ai is the staggered quark mass 
matrix. For n Kogut-Susskind (KS) flavors, is a Anrir x Anrir matrix. 

• \ 



M 



( muuI®lR rnudI®lR rriusl ® Ir 

rrid J ® Ir niddl ® Ir nidsl Ir 

rrisul ® Ir rUgdl ® Ir rugsl ® Ir 

V ; ; ; 



(B.3) 



/ 



where 7 is a 4 x 4 unit matrix, and Ir is rir x unit rephca matrix. In short, 

M. = m.® I ® Ir, where 



/ rriuu ruud rrius 
rridu rudd ruds 



m 



rrisd iris 



\ 



\ 



J 



(B.4) 



is the n X n quark mass matrix. We know from QCD that the quark current 
dr{x)ur{x) is given by 

l^yM - -^^y (B.5) 

where >Cqcd is the QCD Lagrangian. 

Since staggered chiral lagrangian (jCsxPt) is an effective equivalent Lagrangian 
for >Cqcd in low energy limit, the effective current Ur{x)ur{x) is obtained from 



dr{x)Ur{x) 



dMdrUrixy 



where 



(B.6) 



(B.7) 



is the staggered chiral Lagrangian [60] . We omit the high order terms and the terms 
that are independent of A^, /tt is the tree-level pion decay constant (/jr = 131 
MeV [60]), 5o is the constant with the dimension of the mass [94], and E = 
exp 



[60]. The held $ = X] o ® ^ described in terms of the mass 



6=1 



eigenstate field (f)^ [60] , where 0^ is a 3 x 3 pseudoscalar matrix with fiavor compo- 
nents with fiavor /, /', the index of the taste replica r, r', and taste h which 
is given by generators = {^b,^^bi^Ji,^i^Jiiv,i^,,iI} [60]. Hence, * is 4nn^ x AnUr 
pseudoscalar matrix in SxPT [60], and the subscripts u, d denote its valance fiavor 
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component. Therefore, S is also a Anrir x AnUr matrix. The Tr is the full AnUr x AnUr 
trace. Therefore, the effective current is [94] 

dr{x)Ur{x) = BoTTtMxfU^Ur, (B.8) 

where the notation TVt stands for the trace over taste. 

B.2 Bubble Contribution 

In this section we compute the bubble contribution from two intermediate states 
to the oo correlator in Eq. (B.l). Prom the discussion in Sec. B.l, the point scalar 
current can be described in terms of the pseudoscalar field $ by using SxPT [94] [63] , 

dr{x)Ur{x) = SoTVt[$(x)2]d.„. (B.9) 
Ur{x)dr{x) = BoTT^mx)X^dr, (B-IO) 

where Bq — — ^ is the coupling of the point scalar current to the pseudoscalar 
field which can be numerically determined from the lattice simulation. 

Now we consider the lattice simulations with staggered KS valence quarks and 
three flavors of staggered KS dynamical sea quarks with rriu = rrid ^ mg. We 
only consider the correlators where the valence-quark mass is equal to the u or d 
dynamical sea-quark mass. We know that our oq correlator is estimated by using 
the taste-singlet source and taste-singlet sink. This ao correlator takes contribution 
from the ao meson and the bubble contribution. For concreteness, the bubble 
contribution to the ao correlator in a theory with rir tastes per flavor and three 
flavors (namely, iVj = 2 -|- 1) of KS dynamical sea quarks is [60] 

K^'^^'i^) - {dr{x)Ur{x) uAO)dAO))^^\ 

r \ ) 

= — E l{'^t[^\rUrTr,[^\,d^,)], (B.ll) 
Ur r,r'=l { ) 

where the subscript ao specifies the bubble contribution for ao meson. If we consider 
this identity TVt (t^'T^) = ASab, we arrive at 

tj2 16 Nf rir 

Bif\x)^-^Y. E E {<t^iu{^)<t>iuM <Pij,mu,m. (B.12) 

"'r o=l i=l r,r' = l 
6=1 J = l t,t'=l 
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where a. 6, a', h' are taste indices, i.j are flavor indices, and r, t, r' , t' are the indices 
of the taste replica. The Wick contractions result in the products of two propagators 
for the pseudoscalar fields. Therefore, the bubble contribution can be expressed in 
terms of the pseudoscalar propagators {(f)^ (ffi) . 

jd2 16 rir rir ( 



where we use charge symmetry to equate correlators with u d. 

First of all, we sum over the indices of the taste rephca in Eq. (B.13). For 
illustration, we consider the second term in Eq. (B.13) for taste b — I,V, A. 

r,r' t,t' I ) 

- e|««.(^x,..,(o))(<^i..(^x,u.,(o))| + 

r,r' \ ) r^r' 

r,t I J r^t 

r,t I J r^t 

The first term and the second term in Eq. (B.14) should be zero, and the third 
term in Eq. (B.14) should be 

Uririr - l)(0t«. <^t«.)cc (^td. </'l«.>cc, (B.15) 

where the notation CC stands for the connected contribution. The fourth term in 
Eq. (B.14) should be 
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{<t>luS^) <t>luX^)hc{€.aX^) 0l«.(O))cc|, (B.16) 
where the notation DC stands for the disconnected contribution. Hence, the second 
term in Eq. (B.14) can be changed into 

r,r' t,t' I ) 

nU€.uM €ruM)cc (0))cc + 

0l«.(O))Dc(0l,.(x) 0l„,(O))cc. (B.17) 

Because there is no disconnected contribution for taste b ^ /, V, A, therefore 
Eq. (B.14) can be written as 

r,r' t,t' { ) 

= nl{ct>luM </'U(0))cc(C(x) 0l„.(O))cc. (B.18) 
Prom the similar procedure, for taste 6 = /, V, A, we obtain 

EE(«d.(^) A,u,Xm€.uM) <,..,(o))| 

r,r' t,t' I J 

= nr{<Pl,^,{x) 0l„,(O))cc(0U(x) <,,(0))dc, (B.19) 

where we consider the 0^^) does not include the connected contribution. For 
taste 6 7^ /, y. A, this term should be zero. For any taste 6, we obtain 

ee(«.*(^) C(o))(C(^) <,s.,(o))| 

r,r' t,t' I J 

= </'U(O))cc(0U(^) C(0))cc. (B.20) 

From now on, since there is no summation over the index of the taste replica. 

We suppress it. Thus we can rewrite Eq. (B.13) as 

16 r 



6=1 I 

{<f>li^) C(0))cc (C(^) 0L(O))cc} 

+ 45o' E 0L(O))dc 0L(O))cc|, (B.21) 

The propagators (0'' 4>'^) for the pseudoscalar field 0^ for various tastes 6 are studied 
in Ref. [60]. Propagators for all tastes (/, V, A, T, P) have connected contributions. 
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while only tastes /, V and A have disconnected contributions. Hence, Eq. (B.21) 
can be rewritten by plugging in the mesonic propagators, that is, 

B!f^{t) = F.T.[B!f''{p)]^=o with (B.22) 
BS.pt. . ^ Bl i 1 k^ + Ml 



16 



11 11 

(A; + p)2 + A;2 + M^^ ^ {k + pf + M^^ A:^ + M|.^ 



6=1 

1 12 1 1 



6=1 

~^ (A; + p)2 + A;2 + M^^ + 3 (A; + p)2 + M^^ A;2 + M^^ i"' 



where 5y = o^f^y, 5a — a^f^^ [63] [60] 



1 . 2 



M,, = -M2^ + -M|^ (B.24) 



3 3 

m: ' 



M. 

Here 



.V = 2G^^v + ^i + 3My-^y) (B.25) 
?y - i(M^^ + M|^+3My + Zv). (B.26) 



= ^(MJ^ - M2 J2 _ 2nr5vm^ - M^^) + 9My)2. (B.27) 

For the axial taste we just requires V ^ A. The pseudoscalar masses Mu^ = M^^ = 
2Bomu + a2 Ab, = M^s, = 2Soms + a'^Ab and M^, = M^^^ = 5o(m„ + m^) + a^Af, 
have been determined by MILC simulations [38, 10]. The taste breaking term a^Ab 
is independent of the flavor [60] and will disappear in the continuum limit (i.e., 
a — )■ 0) . We can prove that 

k'+Ml 1 f 2 Cy^ Cy^, \ 

{k''+M^^){k^+M^y){k-'+M^,y) ArirSy [k^+M^^ k^+M^y k^+M^>y J ' 

(B.28) 

where 

(Ml -M?r -nJ.. + Z,A 
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^v„, — 



{Ml - - nM - Zv) 



(B.30) 



For the axial taste we just requires V ^ A. We simplify Eq. (B.23) as 



-8 



+ 



4Ci 



(A; +p)2 + k^ + {k+pY + M^^ k^ + M^y 



+ 



4Ci 



Vr, 



(k+pY + M^^k' + M^^ {k+pY + M^^k^ + M^^ 



16 



h- 



1 



+ 



1 



U (k+pY + M^^k' + M^^ {k+pY + Ml,^k' + MlJ 
1 12 1 1 1 



2 

+ 77- 



{k+pY + M^^k^ + M^^ ^{k+pY + M^^k^ + M^iY 
For 1 + 1 + 1 theory, we obtain 



(B.31) 



T3^ 1 Q 



2e 



-(^(p+k)2+M2^ + ^k2+MiJ^)t - (^(p+k)2+M2^ + ^k2+M2^) t 



3 - /^p+k)'+M2 Jk2+M2, 



- (^(p+k)^+M2^+Vk2+M^^) t 



p+k)^+M2^^k2+M2^ 

(^(p+k)^+M^^+^kHM^V)* 



+4Ca, 



(p+k)^+M2^y'k2+M2^ 

;p+k)2+M^,^+^k2+M^2^)t 



^(p+k)^+M2^ykH]i^ 



-4C 



(p+k)^+M^.^ + ^k2+M2,^jt 



p+k)VM^^^/k^TM^ ^(p+k)VM^^yk^TMj; 

(^(p+k)2+M2^ + ^k2+M2^)t -(^(p+k)2 



^^+M^^ + ^k2+M2^)^ 



^(p+k)VM^^^k2+M2^ ^(p+k)VM2^^k^ 



2 



1 16 



-(^(p+k)VM2^ + ^k2+M2jt _^ ^ ^- (^(p+k)2+M|,^ + ^+M^, )t 



8^=1 V(p+k)^+M2^^kHMB; 16 ,=1 V(p+k)^+Mi^ykHMj; 



. (B.32) 



APPENDIX C 



BUBBLE CONTRIBUTION FOR /o 
CORRELATOR 

The mass of the /o meson can be rehably determined on the lattice. In order 
to determine the mass of the /o meson, we use lattice simulations to evaluate the 
scalar /o correlator in Eq. (8.3), that is, 

C{t) = ^(M,(x,t)M,(x,t)M,(O,0)M,(O,0))cC + 

X 

2nr ^{Ur{x, t)Ur{x, t)Ur{0, 0)Ur{0, 0))dC, (C.l) 

X 

where r and r' is the indices of the taste replica. The extraction of the mass of 
the /o meson (J^ = 0^ and / = 0) is straightforward. However there are many 
multihadron states with = 0^ and 1 = can propagate between the source 
and sink of scalar /o correlator. Of special interest is the bubble contribution B, 
which gives a considerable contribution to the scalar /o correlator, and it should be 
included in the fit of the lattice correlator in Eq. (C.l), 

C{t) = Ae-'^fo' + B{t), (C.2) 

where we omit the unimportant contributions from the excited scalar /o meson and 
other high order multihadron intermediate states. 

C.l Bubble contribution 

In this section we compute the bubble contribution to the scalar /q correlator in 
Eq. (C.l). From the discussion in Sec. B.l, the point scalar current can be described 
in terms of the pseudoscalar field $ by using SxPT [94] [63] 

Ur{x)ur{x) = ^oTrt [$(a;)^]„^„^ , dr{x)dr{x) = 5oTrt[$(x)^]d^rf^, (C.3) 

where Eq is the coupling of the point scalar current to the pseudoscalar field ^{x). 
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We know that our scalar /o correlator is estimated by using the taste-singlet 
source and taste-singlet sink. It takes contribution from the scalar /o meson and 
the bubble contribution. For concreteness, the bubble contribution in a theory with 
rir tastes per flavor and three flavors of KS dynamical sea quarks is [60] 

o2 nr c "I 

(C.4) 

where $ is 12??,^ x 12??,^ pseudoscalar matrix in S^PT [60], the subscripts u, d denote 
its valance flavor component, and the subscript /o specifies the bubble contribution 
for the /o meson. If we consider this identity Trt {t°'T^^ = A6abi we obtain 

T^2 16 16 rir rir 

"'r a=lb=li=lj=lr,r'=lt,t'=l 

TD2 16 16 Ur Ur 

ttEEEE E E(C.(^)C...(^)0U(o)4c^.(o))-(C-5) 

"•r a=l 6=1 i=l j=l r,r'=l t,t'=l 



where a,b,a',b' are taste indices, and i,j are flavor indices, and r,t,r',t' are the 
indices of the taste replica. Prom similar procedure in Sec. B.2, we obtain 

16 



1 

b=l I 

i't'Kx) 4'''„m {4'Lix) K.(o)>} + 



''^r b=I,V,A I 



4n,(0l(x) 0L(O))dc '/•L(0))cc|, (C.6) 

where we consider the u, d quarks are degenerate in mass. The bubble contribution 
in Eq. (C.4) can be rewritten by plugging in the mesonic propagators, 

Bff^it) = F.T.[ Bff^ip) ]p=o with (C.7) 
o5xPT/ x_ :B|^r 1 1 1 



{k + pf + M^y k^ + M^y ^ (A; + pf + M^,y k^ + M^iy 
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fi2 



(A; + p)2 + M^^ A;2 + M^^ ^^^^ (A; + pf + M^, ^ /e^ + ^ 



-4 



1 



1 



+ 



{k + py + M^y + M^,y {k+py+M^,yk^ + M^y_ 

11 11 

[{k+pf + M^^k'^ + M^,^ ^ Jk+pfTM^T^TM^ 
' -4- ' 



(A; +p)2 + F + M^^ {k+pf + M^^ F + M^^ 
11 11 



16 



6=1 



+ 



+p)2 + k' + {k+pf + Ml^ F + MiJ j ' 



(C.8) 



where M^,, M^^, M^,^, 5^, M^^, M^,^, C^,^, C^^,, C^,, Ca^, are denoted in 
Appendix B. The Fourier Transform of Eq. (C.8) (1+1+1+1 theory) is 



+c\ 



^7 



^2 . 



-CAr,CA, 



;p+k)^+M2^^kHM2^ ' ^(p+k)^+M2^^k2+M2,, 

^(p+k)2+M2^+Vl?TM^)t ^ ^-(^(p+k)2+M2,^+^k2+M2,^)t 

p+k)^+M2^^k2+M2^ ^ ^"'Tl^+kf+^IyPT^ 

- (^(p+k)2+M2^ + ^k2+M2,^)t - (^(p+k)2+M2,^+ VkHM;^^) t 



;p+k)^+M2^Jk2+M2^ 



+ 



;p+k)^+M2^Jk2+M2^ 



^-(^(p+k)^+M^^,+^k^+M^^,,)t ^- 



(^(p+k)2+M2,^+ 



+ 



;p+k)^+M2^Jk2+M2^ ^(p+k)^+M2^Jk2+M2^ 



-4- 



-(^(p+k)2+M2^+^k2+M2^)t -(^(p+kr 



^(p+k)VM2^^k2+M, 



-4- 



2 



(p+k)VM2^^kHM, 



2 



1 

-E 



-(^(p+k)2+M2,^ + ^/k2+M^^ It 



-(y'(p+k)2+Af2^ + y'kHM2j, 

^(p+k)VM2^y'k2+M2^ ' ^(p+k)VMi^^k2+M: 



+ 



(C.9) 



where Bl = ° 
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APPENDIX D 



LBO TREATMENT OF DECAY CHANNEL 



FOR HYBRID 



In this dissertation, we are interested in the decay of a hybrid exotic quarkonium 
l""*" state H to Xb state with the emission of a pair of pseudoscalar mesons in a scalar 
channel (tttt). To make the problem simple we assume that the scalar meson (tttt) 
in the final state does not interact with the Xb final state. The lattice calculation 
is performed by the LBO approximation, and the wave functions of the initial and 
final state are easily described in the rest frame of the heavy quarks. We call it 
the "body- fixed frame." The displacement vector R of the heavy quarks specifies 
the z axis of the body-fixed frame. The exotic hybrid state (H) is n„ and the final 
quarkonium state is E+ state (x&). 

As discussed in Chapter 4, the potentials are computed in the static approxima- 
tion. The wave function of the hybrid exotic in lab frame is ^ 



where we ignore the spin degrees of freedom, m is the projection of the total angular 
momentum (excluding spin) on the lab z axis, and ±1 is its projection on the 
body-fixed axis. The Euler angle a specifies a rotation about the body-fixed axis 
and it is an internal coordinate for the excited flux tube. 

^Z/^-normalized Wigner D-functions Dj^ f^,{(j),6,a) and ordinary Wigner D-functions 



•^h{R. 0, e, a) = ph{R) 



m,—l 



{(f), 9, a) 



(D.l) 



-^M M'{4'j ^1 ct) have relationship, D'lj 




/ da sin 6d6 



Jo Jo Jo 
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In lab frame the state H can be described by 

\H) = j d^-Rd'v I da^^^I^HiR, <t>, e, «) Cj(ri)C^(r2) |0) , (D.2) 
where ri = (R - r) /2, ra = (R + r) /2, and 




(D.3) 



here the o^a is the creation operator of gluon. 

By convention the wave function of the Xb state in rest frame is 

^^,{R,cl>,e)^p^,{R)Y^^,{e,cl>). (D.4) 

In lab frame the state A can be described by 

\A) = / dWr^^*^,(i?, 0, 0) 4(ri)4(r2) |0) . (D.5) 

In the body-fixed frame with R fixed the initial hybrid exotic state has only the 
internal gluonic degree of freedom with a wave function cosa/y^ corresponding 
to the choice A = ±1. We can model this transition in the lab frame through an 
interaction potential 

Vt" = |d3Rci3r|rf='r;x(R0Cj(ri)4(r2)CQ(ri)CQ(r2)x 

V'V'(rs + r) [wi,i{r'Jag,i + w/,_i(r;)ag,_i] , (D.6) 

where the primed coordinates arc the coordinates of the scalar meson in the body- 
fixed frame, wj is the wave function of the scalar meson in the body-fixed frame 
as determined by the dynamics of the transition, and x{R) controls the transition 
amplitude. Since final state Xb has A = 0, the scalar meson must carry off the 
nonzero ^-component of the angular momentum in the body-fixed frame. 

The wave function of the scalar meson is a free spherical wave with an appro- 
priate angular momentum {L — 1). Hence, its wave function in rest frame is 

i/jsivs, Os, (t>s) = Ns3i{pr)Yi,mX^s, (ps), (D.7) 

where Ng is a normalization factor, ji{pr) is the spherical Bessel function, and 
^i,ms(^s)0s) is the spherical harmonic function, and p is the magnitude of the 
angular momentum. 
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In this notation we evaluate the transition amphtude in lab frame by 

cos(q; — 0^) 



271 



X 



7r(pi)7r(p2) ^V'(O) O) S^^\Ph -Pa- p), (Di 



where |7r(pi)7r(p2)) = flpi^pa |0) ( namely, flp^flpj is the creation operator for tttt), 
p = Pi + P2, and the p is the magnitude of the angular momentum p. For brevity, 
we divide (^4, 7r(pi)7r(p2) V}^^ i?) into two parts, that is, 

(D.9) 



A7r(pi)7r(p2) V;-" h)=x^^{p^,P2)6^^\Ph-Pa-p), 



where 



X ^^^^^=^x(R)^,«) (7r(pi)7r(p2) l^l^m 



. 



(D.IO) 



The integration over scalar meson coordinates is carried out in body-fixed frame. 
To this goal we use the rotation property of the spherical harmonic function 



Y,,m. {0s, 0.) = E Dln^A^, 0)Y^M, 
A 



(D.ll) 



Because the volume element is rotationally invariant, we have 

(pi,P2) = j RHRdn j r'^dr'sdn'^Y.D]:,,M9)x{R)p^SR)pH{R)N,^^^^^^ 
X yi,A(^;, ct>',)Y*^^{9, <l>)x{R)wj(r's, 9',, 0',) (7r(pi)7r(p2) ^^^'(0) 
X j da 



cos(a - 1^^,1(0, 9, a) + £>^,_i(0, 9, a) 

a/2 



0; 

(D.12) 



After integrating out the a and rearranging terms, we obtain ^ 

xUPi,P2) = j R^dRdn j r'^dr'^dn'^Y.D]:^A^4>,0)p^SR)pH{R)Ns3i{pr', 

X yi,A(^;, 0'jy4,(0, 9)x{r)wj{t',, ^'j (7r(pi)7r(p2) m^) 



^Here we use the convention that Z)4 n{4'i ^1 Qi) = e'™'^d^ ^, a)e' 
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(D.13) 

When carrying out the 0^ integration and rearranging terms, we can omit the 
summation over A. In order to normahze correctly later, we here still keep the 
integral over 0^, that is, 

?Vx.(^)a^(P,^)PH(i?)rfi?x (7r(pi)7r(p2) |^V^(0)|0 



where 

x{p,R)=x{R) I rf(ir^cos^^d^^Ar,ji(pr^)yn(^^,</''J^/;7(r^,^^>'J. 
We can divide Eq. (D.14) into three parts, that is, 

XnAPi^P2)^ j R'^PH{R)x{p,R)py^^{R)dR X F{m, 171^,7111,), 

where 



(D.15) 



(D.16) 



F(m,m„m,) = J dnDl:^^,(<l>,e)Y*^^(<l),e)e^^^[dm,m - (D.17) 

We choose the right measurement choice of spherical harmonics function Y and 
Wigner D function, that is. 



and 





V Stt 




(D.18) 


YM'l',0) = 


a 

\ / — cos 
V 47r 




(D.19) 




V Stt 




(D.20) 






1 -cos^ 
2 


(D.21) 






_ l + cos^^^,^ 
2 


(D.22) 






sin 9 


(D.23) 


only when 








m = rus 


+ rub, 




(D.24) 



F(m, 777.5,7716) is nonzero. Here we list all the cases, that is. 
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Case 1 
Case 2 
Case 3 
Case 4 
Case 5 
Case 6 
Case 7 



m = 1, rris = 1, nib = 0: F{m,ms,mb) = | 
m — 1, rus — 0, ruf, — 1: F{m, rus, rrib) = | 
m — 0, rus — 1, rrif, — —1: F(m, m^, rrib) — \ 
m — 0, nis — 0, nib — 0: F{m, m^, m;,) = 
m — 0, nis — —1, = 1: F(m, m^, m^) = | 



m — —1, = —1, mft = 0: F(m, m^, m^) 
m = —1, nis = 0, mfe = —1: F{m,ms,mb) 
Finally, we obtain 

xUpi,P2) R^Ph{R)x{p, R)p^,{R)dR x (7r(pi)7r(p2) |V^V'(0)| o) . (D.25) 

If we introduce the reduced radial wave function uh{R) = Rph{R), and the reduced 
radial wave function u^^{R) = Rpy^^{R), we rewrite above equation as 

'o). (D.26) 
(D.27) 
(D.28) 



Xt,^{Pi,P2) = 2 / uh(R)x(p, R)u^f,(R)dR x (7r(pi)7r(p2) V'V'(O) 
If we denote 







x{p) = ^/ uh{R)x{p, R)u^^,{R)dR, 

we have 

x„n{p) = x{pi,p2) (vr(pi)7r(p2) ipipiO) 

In Chapter 8, we use Eq. (D.26) to calculate the Xt^t^^p). 

The integration is done in the body-fixed frame. Since the transition amplitude 
is computed on the lattice, we replace the continuous integration with a lattice sum. 
The wave function is measured in lattice momentum space for a lattice of spatial 
size a^L^ through the expansion 

(D.29) 



wi{x,y,z)^ ^ wi{kx,ky,kz)sm{kxx)cos{kyy)cos{kzz), 

In this approximation 

^ Yl ^^^^j^(pK)yn{0's,(l>'s)x(k,R) sm{k^x)cos{kyy)cos{k^z). 



x{p, R) 



L3 



(D.30) 

With a sufficiently large lattice volume the overlap sum over coordinates {x, y, z) 
recovers the continuum momentum conservation constraint p = |k|. Here we use 
Eq. (7.22). 
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Now we discuss how to calculate the above equation numerically. First we 
calculate the Ng, this number come from the normalized condition 



x,y,z 



where 



Yu{9,cf>) = 



X + iy 



, , sinOpr) cos(pr) 

:J^(p^) = 



(D.31) 

(D.32) 
(D.33) 



p^r^ pr 

here r = y/x^ + + z^- Then we can use Eq. (D.30) to calculate x{p, R) at any 
given momentum p. 

In practice, we can rewrite Eq. (D.25) as 

x..{puP2) = E E a'N,Mpr',)Yn{9:, 



where 



X x(k) (7r(pi)7r(p2) i>i>{0) O) sm(kxx) cos{kyy) cos{kzz),(D.3A) 



x(k) = ^ / UH{R)x{\i,R)u^^{R)dR. 



From Chapter 7, we know, 

V(pi)7r(p2) V^(0) 



0) = 6. 



^2K(pi)V2^.(P2)' 
where ^tttt is denoted in Eq. (7.42). For notational simplicity, we define 



(D.35) 



(D.36) 



(D.37) 



If we define 



^TTTrb) = ^ E E (t>'s) ^(k) s\n{kxx) cos{kyy) cos{k,z), 



then 



^2K(Pi)V2K(P2)' 



(D.38) 
(D.39) 



APPENDIX E 



THREE-BODY PHASE SPACE 

In this dissertation, we are interested in the decay 

^^Xb + TT + TT, (E.l) 

where two-pion final state is in an scalar singlet channel. In Appendix D, we 
calculated the transition amplitude x^7r(p) for two-pion production as a function of 
the momentum p of the two-pion system in the rest frame of the H. 
The starting point is the Fermi's Golden rule, that is, 

dV = j^^\{f\V,\i)\\2nY5\pf-p,)df, (E.2) 

where {f\Vi\i) is the transition amplitude from initial state i to the final state /, 
df is the phase space factor. 

The general rule is that the various factors in the phase space are determined 
by the normalization of the states. Prom the definition in Eq. (D.2), we can show 
that the initial state \H) is normalized as 

/rf='pH|pF)(p^f| = 1, (E.3) 
and from the definition in Eq. (D.5), the final state \A) is normahzed as 

j d''pA\pA){pA\ = \ (E.4) 

Hence, they should not give a factor in the phase space formula. There are two 
other momentum factors in our phase space formula, called pi and p2- From the 
Ref. [64], they are normalized as 

|PiP2 ) ( P1P2I = 1, (E.5) 

pi P2 
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They are replaced hjp = pi2, the total momentum of the two pions and pi = k, the 
momentum of one of the pions. In the analysis of the bubble diagram the momenta 
p and k are quantized in a box and summed. Hence the phase space factor 

where £'i,i?2 are the energies of two pions. Hence, according to Fermi's Golden 
rule, we obtain transition rate (F) by 

S(\Pi2\-p)x^Ap)^, (E.7) 

where pn is the four momentum of initial hybrid exotic state H, pA is the four 
momentum of Xh final state, pi is the four momentum of one pion final state, p2 is 
the four momentum of other pion final state, we denote pi2 = pi +P2, and where 
x-ktt{p) is denoted in Eq. (7.40), which is the transition amplitude from initial state 
to the final state. 

To simplify the integration in Eq. (E.7), we introduce the integral over the 
intermediate four momentum pi2 and its invariant mass s = p\2 — £■12 ~ |Pi2|^- We 
rewrite the Eq. ( E.7) as 

r = (^(^j J dpdEi2d^PAd'pi2{2n)^S\pA + Pi+P2-PH) 

X Fi2(s)5(|pi2| -p)£..(p)', (E.8) 

where d^pA — |paPc?|pa| smOdOdcf), and ^12(5) is the two-body phase space factor 
for two-pion system, that is, 

i^l2(«)=/^|f<^V2-Pl-P2), (E.9) 

where d^pi — |pipd|pi| sin This phase space factor is a Lorentz invariant. 

Thus we can easily evaluate it in the rest frame of two-pion system, that is, 

I 1 2 

F,2{s) = ^7: j ,^^d\p,\5{E, + E,-^s) 
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where m^r is the mass of the pion meson. Therefore, we obtain 

r = 27r' i^—j J dp dEi2ct^pAd^pi2S\pA +P1+P2- Ph) 
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In the rest frame of the H we integrate above equation, and we obtain 



r = 8.^ C^S 



\27rJ 



// 1 771 
dpdEi2\pA\^d\pA\S{EA+Ei2-MH) J 2' ^ ^(\Pa\-p)^^-^p)^' 

(E.12) 

In the rest frame of the H we integrate over |p^| — p and s to get 



dpp x^^(py 



mt 



\A {AM^-p^' 



(E.13) 



where we approximate Ea — Ma, and define AM = Mh — Ma- If we consider 
(AM)^ -p^ > Ami, we have 



p^ax = ^(AM)2 - Ami. 



(E.14) 
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